A FUZZY SETS APPROACH TO INVENTORY 

SYSTEMS DESIGN 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the degree of 

MASTER OF TECHNOLOGY 


By 

RAJESH DVBE 


to the 

INDUSTRIAL AND MANAGEMENT ENGINEERING PROGRAMME 
INDIAN INSTITUTE OF TECHNOLOGY. KANPUR 

APRIL. 1977 



ii 



gertifiqatb 


4 ? 


Shis is to certify that this work entitled 
Fuzzy Sets Approach to Inventory Systems Design^ 
has been carried out under my supervision and that 
it has not been submitted elsewhere for the award 
of a degree . 


April, 1977 



( J.L, Batra) 


Convener \ 

Industrial and Management EnggL Programme 
Indian Institute of Technology, 

Kanpur 208016 


l 

I ! " ■■ , fi v--’ r 

\ ^ ^ I 


Vv ’ 


!.,,\vAn 


. -.. a. I- ^ 



f FUR 

Cl' I . ! jl:iU I RY 

smg 


1 1 AUG I37f 




TMEP- Dufe-Fuz 



iii 


AGI2T0WLEDGEMEITTS 

The acknowledgement usually . tiirns out to be the 
most pleasent part of a thesis to write- This is not only 
because it is the last page to be written, but primarily 
because it is the place where the author can make explicit 
his debt to the many people who have helped to make the woik 
possible, I wish to thank; 

- Dr, J,l. Batra for initiating me into this topic 
and for his constant encouragement, help and 
guidance. 

- Dr. R, Raj agopalan whose work in this area has 
been a constant source of inspiration to me. 

- Dr. J.K. Dwivedi and Mr, S.R. Shrivastava of the 
Health Center, IIT Kanpur for help in data collection. 

- Drs. J. Bezdek, J. Dunn, R. Bellman, I. Gitman, 

R. Gapocelli, A. De luca and E. Ruspini for 
promptly acceding to my requests to send their 
papers in the field. 

- Messers H. Karnick, V. Ghanohlani, Priya Ranjan 
and V, Purohit who helped me in several ways. 

- Mr. J.K. Misra for immaculate typing, and 

Mr. Buddhi Ram Kandiyal and Mr. R.S, Tripathi for 
cy do styling. 


Rajesh Dube 



IV 


OONTENTS 


Chapter Page 

LIST OP TiiBLES vi 

LIST OP PIGURES vii 

SYNOPSIS viii 

I. INTRO LUGTION 1 

1.1 The Item Glassification Scheme 5 

1.2 Porecasting the Requirements 6 

1.3 Lot-Sizing 7 

1.4 Reorder Levels 8 

1.5 Organisation of the Thesis 9 

II. LITERilTURB REVIEW 10 

2.1 ABC Analysis 10 

2.2 Puzzy Sets 12 

2.3 Cluster Analysis 14 

2.3.1 Hard Clustering Techniques 14 

2.3.2 Puzzy Clustering Techniques 15 

2.4 Porecasting 17 

2.4.1 Quantitative Porecasting Methods 17 

2.4.2 Qualitative Porecasting Methods 19 

2.5 Inventory Control Techniques 20 

2.5.1 Analytical Techniques 21 

2.5.2 Graphical Techniques 27 

2.5.3 Lynamic Models 28 

2.5.4 Simulation Models 30 

2.6 Computerised Inventory Systems 31 

2.7 Summary 35 

III. PROBLEL'^I STATEMENT AJD METHODOLOGY 37 

3.1 Statement of the Prohlem 37 

3.2 Inventoiy Management Objectives 38 

3.3 The Hospital Medicine Inventory 

Control Problem 39 

3.4 Selective Inventory Management Concept 41 



3.4.1 Traditional Analysis 42 

3.4.2 Tlie ¥eed for a New Glassification 

Methodology 43 

3.5 The Fuzzy Sets Approach 44 

3.5.1 Fuzzy Sets and Fuzzy Clusters 46 

3.5.2 Importance Through Feature Vectors 55 

3.5.3 Glassification by Fuzzy ISODATA Process 65 

3.5.4 Ghoo sing the Initial Centroids 70 

3.5.5 Measures of Partition Quality 71 

3.6 Selecting the Ordering System 78 

3.7 An Overview of the System 80 

3.8 Control Procedure for 'Class I' 

Medicines 83 

3_.8.1 The Forecasting Scheme 83 

3-.8,2 Lot-Sizing Procedure 86 

3.8,3 Calculation of Reorder Point 88 

3.9 Control Procedure for 'Glass II' 

Medicines 93 

3.10 Control Procedure for 'Class III' 

Medicines 94 

3.11 Summary 96 

A GilSE STUDIs RESULTS, DISCUSSIONS AND THE 
SCOPE FOR FURTHER WORK 97 

4.1 ABC Classification 98 

4.2 Fuzzy Inventory Classification 99 

4.2,1 Initial Centroids 99 

4. 2. 2 Partitioning by Fuzzy ISODATA 108 

4.3 ABC Analysis vs. PICS 116 

4.4 Integrated Forecasting and Lot-Sizing 119 

4.5 Computational Experience 119 

4.6 Scope for Further Work 122 

4.7 Summary 125 

ESFEEBNOES 126 

APPENDIX I 136 

APPENDIX II 144 



Vi 


LIST OF TABLES 

TaLle Page 

4.1 Membership Grades of the Medicines in the 

Three Classes, 101 

4.2 Quality of Clusters for Yarious *m' and 

*c^ Values. 109 

4.3 Comparison of ABC Analysis vs. PICS. 116 

4.4 Output of the Computer Program for Computing 

Order Quantity and Reorder Level 120 



vii 


LIST OF FIGURES 


Figure 


Page 


1 . 


2 . 

3. 

4, 


5. 


6 . 


7. 

8 . 

9. 


10 . 


Logic of tlie Fuzzy I SOL AT A Process 69 

Selection of Initial Centroids 72 

General Procedure for Item Classification 77 

General Logic Plow Liagram for Computing 
Lot Sizes and Reorder Levels 82 

Annual Rupee Usage Curve of Sample of 119 
Medicines Stocked at IIT Eanpur Health 
Center 100 


(i) Questionnaire Used for Lata Collection 105 

/■ 

(ii) Form for Lead Time and Lemand Lata 

Collection 106 

Generation of Initial Centroids by Ball 


and Hall Method 

Variation in H and F with Exponent ’m’ 
o c 

for c = 2. 

Variation in H and F with Exponent ’m’ 
for c = 3 

Variation in H and F with Exponent *m' 


107 

110 

111 

112 



for c = 4 



viii 


STUOPSIS 

The present work deals with the development of 
an integrated inventory olassif ication, forecasting and 
control system. This work is the first attempt to cover 
the gap between the concepts of Selective Inventory 
Management and the way it is practiced in the form of annual 
rupee usage based iBG analysis, ilso, it is for the first 
time that the theory of ’Fuzzy Sets' has been used for 
inventory systems design. Using this theory a computerised 
inventory classification .system has been developed. The 
proposed system is quite compatible with the concepts of 
Selective Inventory Management which emphasise the 
allocation of management control effort according to the 
overall relative importance of the items and not simply on 
the basis of annual rupee usage as is done in the classical 
130 classification. 

The overall importance of various items has 
been represented by feature vectors and the ’^uzzy 
ISODATA' algorithm proposed by Bezdek and Dunn has, been 
used for classification. A hospital medicine inventory 
situation has been studied and three groups of medicines 
with varying degree of importance have been found. The 
first group consists of those medicines which are Important 
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fr.om all considerations, the second group has the medicines 
which have low annual rupee usage and high expiry time but 
are otherwise important, and the third group comprises of 
the medicines which are unimportant from all considerations. 

The proposed classification scheme has been found superior 
to the traditional ABO analysis. 

iin adaptive control limit system has been designed 
for the first two classes of important medicines. A forecasting 
technique given by Buchan and Koenigsberg and a modified 
version of Be Matteis^s Part Period Algorithm have been 
used for lot sizing. An extension of Tchebycheff 's inequality 
has been used to set reorder points for specified service 
levels. A more economical static control limit system has 
been proposed for the third group of unimportant medicines. 

A computer program has been developed for all the three 
classes of medicines. 

This approach can be easily extended to cover 
other pure inventory situations like warehouses, grocery shops, 
super markets, etc. and should result into significant 
improvements in the existing cost structures and service 


levels 



CHAPTER I 


IHTROPUCTIOH 

Economic development provides increased output and 
efficiency through, specialisation. In the caveman's society, 
each mah produced all of his own goods - food, clothing 
and a few utensils. In modern industrialised societies, each 
individual produces only a minute fraction of the goods and 
services he requires. Due to the spatial and temporal 
separation created hy this high degree of specialisation, 
such a society cannot operate without large stocks of 
inventories. Inventories provide the buffer which allows 
these separations to exist. They perform the basic function 
of decoupling the operations involved in converting inputs 
into outputs, 

A price has to be paid for everything one gets, and 
inventories are no exception. Inventories tie up assets and 
additional direct charges accumulate because of the storage 
space requirements, insurance, taxes and so on. Pinas, 
therefore, cannot afford to keep unrestricted inventories at 
eveiy point in the production pn>eess despite the be fits 
they bring. Bven with closely controlled inventories, 
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manufacturing fiims often find that as much as 25 to 50 percent 
of their total assets are invested in inventories. Moreover, 
wholesalers and retailers occasionally find as much as 75 to 
80 percent of their assets in inventories. A survey of 1000 
Joint Stock Companies hy the Reserve Bank of India revealed 
that more than Rs. 1000 crores is blocked in inventories 
alone. 

It cannot be overemphasised that decisions regarding 
the inventory are very important. Inventory models and 
mathematical description of inventory systems provide a basis 
for these decisions. Scientific inventory control is concerned 
with using mathematical models to obtain rules for operating 
inventory, systems so as to minimise the sum of the various 
costs involved. 

The present work examines a special type of inventory 
situation. An attempt has been made to devise a practical, 
efficient and readily adaptable system for controlling medicine 
inventory in hospitals. The proposed system can be used with 
some modifications in situations like drug stores and ware- 
houses, supermarkets and retail shops. It is important to 
note that in industry the criteria usually have economic origins 
but in hospitals it is necessary to consider life and death 
as well as economic criteria. Practice most often has been 
to assign economic costs to aH factors. In this presentation, 
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the influences of various factors on maintaining of inventory 
of medicines is considered in a scientific way tliroughi tlie 
proposed model rather than an arbitrarily assigned cost. 

Some of the important problems encountered in 
designing a medicine inventory control system are j large 
number of items, high and varied service levels, considera- 
tion for noneconomic factors and substitutability and 
perishability of items. The situation is further complicated 
by seasonality and trends in the demand pattem, special 
storing methods for some medicines and variability in lead 
time. 

The first phase of the work deals with classification 
of medicines in various importance groups for deciding the 
level of inventory control effort by the management for 
various classes of medicines. The second phase consists of 
the selection of appropriate procedures for forecasting, 
lot sizing and reorder point setting for a (R, Q) type 
system, depending upon the desired service levels and lead 
time variability. 

1.1 The Item Classification Scheme 

¥ith thousands of items in inventory, it is often 
not possible or even desirable to exercise close control 
over all the items. One method frequently used, to determine 
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to decide the amount of attention to be paid to an item on 
this basis. Moreover, the percentage of items to be put in 
a particular class is arbitrary. Highly dissimilar items 
from other considerations for control may be grouped 
together. This reflects upon the need to group items in a 
more general and systematic way. in attempt has been made 
in this work to develop a general approach for classifica- 
tion which takes into account all the pertinent factors 
contributing to the relative importance of the item. 

The vague and subjective nature of the term 'importance' 
has been the prime source of motivation for using the 
recently developed theory of 'Fuzzy Sets* [1] in this work 
for classifying items in various importance groups. The 
fuzzy classification scheme takes into account the impre- 
cision associated with real -life systems. The central idea 
of this theory is the assigning of a membership grade between 
0 and 1 to every member of a fuzzy set. Bach item is 
represented by a feature vector whose components consist 
of various factors to which the importance of the item can 
be attributed. 'Fuzzy ISODATA* clustering algorithm has 
been used to identify the intrinsic structure of the feature 
vector data. The outcome is the three distinct groups of 
items. First group consists of itenjs important from all 
considerations, the second group is of items which have low 
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annual rupee usage but are otherwise important, and the 
third group consists of unimportant items. The results of 
this classification scheme have been compared with those of 
the traditional ABC classification scheme. The information 
obtained about the various groups present in the data 
has been used to develop a procedure for forecasting and 
lot sizing. 

If 2 Poreoasting the Requirements 

A key aspect of any decision making situation is 
being able to predict the circumstances that surround that 
decision and that situation. Such predictions, generally 
handled under the title forecasting, have been identified as 
a key subpart of the decision-making process. A number of 
forecasting methods or techniques have been developed during 
the last two decades. Selection of proper method plays a 
vital role in the successful application of the knowledge in 
this field. 

It is a common misconception that the overriding 
criterion in selecting a forecasting method for a particular 
situation should be the accuracy of the available methods. 
Several other aspects need to be examined before any method 
is finally selected. The characteristics of forecasting 
methods that have been found most helpful in making the 
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judgemental decisions arei time horizon for which the 
forecast is to be made, demand pattern, model type inherent 
in the method, accuracy, applicability, and associated costs 
in terms of computational effort, data collection and 

i 

s to rage re quire me nts . 

The hospital medicine inventory management problem 
requires a method capable of making medium range time horizon 
forecasts for a demand pattern with seasonality and trend. 

The method must be computationally efficient as demands 
for a large number of items are to be regularly forecasted. 

The present work uses a simple method proposed by Buchan 
and Koenigsberg [2j. The monthly forecasts thus obtained 
have been used as an input for determining lot sizes for the 
various items. 

Lot- Sizing 

A large number of mathematical models have been developed 
during the past years to study various inventoiy situations 
but not all of them are of practical significance. Besides 
the classical BOQ model and its variations which use diffe- 
rential calculus approach, simulation and dynamic programming 
have also been employed to solve complex inventory control 
problems. Simulation requires the foimulation of an inventory 
model as a quantitative measure of performance and objectives. 
The relationships existing within and outside the firm are 
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quantified and optimal results are obtained by chsjiging 
those variables over ■which control is possible until the 
best solution is found. lynamic prograxnming can provide 
better results than the static EOQ model when the demand 
is unstable. However, dynamic Programming h^ its own 
limitations in terms of computer time and memory requirements 
and it does not appear to be a good approach when the 
number of items is large. In this work, a modified version 
of the part periods algorithm [3] has been used. This 
algorithm takes care of varying demand and is also computa- 
tionally efficient. However, for unimportant items, the 
classical EOQ formula has been used. 

1.4 Reorder Levels 

The reorder levels are established to cater to the demand 
during lead time. The procedure used in this work takes 
care of the forecast error and variation in lead time. 
Tchebycheff ’ s inequality has been used to detoimino the 
reorder points for desired confidence levels. To be 
economical in exercising the control effort, different 
procedures have been suggested to estimate lead time and 
its variance for different classes of items. 
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1*5 Organisation of the Thesis 

After this brief introduction we shall review in 
Chapter II the literature relevant to the present work. 

' Chapter III is devoted to describing the models contained 
in the present work and the resxilting inventory control 
system. In Chapter lY, the applicability of the approach 
developed in this thesis is demonstrated through a case 
study from the Health Centre, Indian Institute of technology, 
Kanpur. 



CHAPTER II 


LITERATURE REVIEW 


The literature relevant to the present work encompasses 
ABC inventoiy analysis, the theory of fuzzy sets, cluster 
analysis, fuzzy clustering, forecasting, inventory control 
and the design of computerised inventorj^ control systems. 

In this chapter we shall briefly review the important 
developments in these areas, 

2.1 ABC Analysis 

ABC analysis is a basic industrial engineering tool. 
It involves ranking the items of a population in descending 
order of some exhibited activity and then developing the 
most appropriate technique of handling the high-activity 
'A' group of items, with perhaps different techniques for 
handling the medium activity ’B’ items and the low activity 
'C* items. In many cases, such differentiation in 
techniques for physical handling, control, or management 
gives substantially superior results to treating both the 
important and the unimportant items in the same way. 

A B 0 inventory classification was first introduced 
by H.Pord Dickie [4,5*63 of General Electric Company, USA 
in 1951. It is an adaptation of Pareto’ s law. 
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Vilfredo Pareto (1842-1923), an Italian economist and 
sociologist, in a study of the distribution of wealth and 
income in Italy in 1897, observed that a very large percentage 
of the total national income was concentrated in a small 
percentage of population. Believing that this reflects a 
universal principle, he formulated the axiom that the 
significant items in a group normally constitute a small 
portion of the total items in the gmup and tha,t the majority 
of the items are of minor significance. This law has been 
observed to be particularly true in inventory situations. 

In many inventory situations approximately 20 percent 
of the items account for 80 percent of the total cost (unit 
cost times the usage quantity). In the typical ABC classi- 
fication, these are designated as A- items, and the remaining 
80 percent of the items become B’s and C*s, representing the 
middle 30 percent that account for 15 percent of the cost, 
and the bottom 50 percent that account for 5 percent of the 
cost respectively. However, there is no fixed convention as 
to what constitutes the A, B or G category. The proportions 
may vary from situation to situation. A recent paper by 
Herron [7] discusses some industrial engineering applications 
of the ABC-curves drawn from such analysis. 
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There have been hardly any improvements in the rupee 
usage based ABC classification technique since its inception. 
Though the concept emphasises taking into account all 
pertinent factors contributing to the importance of items 
and not to base the analysis only on rupee value, due to 
lack of mathematical tools available for dealing with fuzzy 
notions embedded in terms like 'importance', perhaps no 
improvement could take place. We have for the first time 
in this work tried to develop a generalised approach to the 
ABC type classification problems. Our approach is based on 
the recently developed theory of fuzzy sets and it uses 
fuzzy clustering aJ-gorithms for classification. 

2. 2 Puzzy Sets 

The concept of fuzzy set, first introduced by 
Zadeh [l] in 196 5 » is a generalisation of the concept of a 
conventioneA set. While in conventional set theory, objects 
are allowed only either to belong or not to belong to a set, 
in fuzzy set theory thej'- are allowed to have a 'degree' of 
belongingness ranging from non-belongingness to belongingness 
to a set. A fuzzy set is characterised by a membership 
function which assigns to each object a grade of membership 
ranging between zea:>o and one with respect to the set. The 
notion of fuzzy set provides a natural way of dealing with 
problems in which the absence of sharply defined criterion 
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of class membership is the source of imprecision which 
makes them immune to conventional non-fuzzy techniques of 
solution. Fuzzy sets, according to Zadeh [8], shall prove 
to have valuable applications in management science, 
economics, information retrieval, artificial intelligence, 
medicine, biology, linguistics and psychology. 

The past decade has seen ra.pid developments in the 
theory of fuzzy sets. A fanrly exhaustive introduction to 
the theory and applications of fuzzy sets [9] lists 238 
papers dealing with a great variety of recent investigations. 
Bellman and Zadeh [10] have illustrated the applicability 
of fuzzy sets in decision making with fuzzy goals and cons- 
traints. Gapocelli and Deluca [ll] have examined fuzzy 
sets in the context of decision theory. Goguen [12] wa,s 
the first to introduce ,'1-Fuzzy Sets* and the concept of 
'entropy’ of fuzzy sets was first defined by De Luca and 
Termini [13] • In reference [14] De Luca and Termini extended 
the concept of entropy to L-fuzzy sets. Another paper by 
De LuCa and Termini [15] discusses many algebric properties 
of fuzzy sets. Rosenfield [16] and Yeh et al. [17] deal 
with fuzzy graphs in clustering. Other useful references 
on the subject include books by Kaufmanh [18] and 
Zadeh et al. [9] and papers by Zadeh [19> 20], Bezdek [21,22,23], 
De Luca and Termini [24] and Gapocelli and De Luca [25]. 
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2*3 Cluster Analysis 

There have been several applications of fuzzy sets in 
devising clustering e.lgo rithms . The primary objective of 
clustering techniques is to partition a given set Of data 
into homogeneous clusters such that all points in the same 
group are similar to each other and dissimilar to those 
in the other groups. With the emergence of fuzzy algorithms, 
now all clustering techniques can be divided into two main 
classes, namely, hard (i.e. conventional) clustering methods 
and fuzzy clustering methods. Hard clustering techniques assign 
each object to one and only one of the clusters, but fuzzy 
clustering techniques assign a membership grade between 0 
and 1 in every cluster to each object. Higher the member- 
ship grade of an object in a fuzzy cluster, stronger is its 
claim for belonging to that cluster. 

2.3.1 Hard Clustering Techniques 

The partitions generated by clustering, irrespective 
of whether they are ’hard’ or ’soft’, can be used in two 
ways. Either the partitions m^ be used to investigate 
the existence of families as done in numerical taxonon^r, 
or they may be used as categories for classifying future 
data points as done in pattern classification. 
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Some of the important contributions dealing with 
taxonomical a,pplications of clustering are by Sokal and 
Sneath [26], Jardine and Sibson [27], and Clifford and 
Stephenson [28j. A book by Duda and Hart [29] gives a 
detailed account of clustering techniques relevant to 
pattern classification. However, there are not many books 
dealing exclusively with cluster analysis. Two of these, 
by Anderberg [50] and Hartigan [51] » ai'e application 
oriented and give a detailed account of computer based 
clustering techniques. An extensive bibliography has been 
provided by Duran and Odell [52]. A paper by Kennedy [55] 
deals with applications of cluster analysis in industrial 
engineering. 

2.5.2 Huzzy Clustering Techniques 

Huzzy clustering is a relatively new area. The 
motivation for adopting fuzziness in clustering comes from 
the fact that in real life, the data is rarely so distinct 
that every member of the data set coiild be described as a 
full member of a single class. A fuzzy model enables us to 
overcome this objection by allowing every individual to 
have partial membership in all the clusters. 

Clustering with fuzzy sets was first proposed by 
Bellman, Kalaba and Zadeh [54]. The papers by Chang [55], 
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Wee [36] and G-itman and levine [37] discuss some of the 
earliest fuzzy pattern classification schemes. In 1969 
Ruspini [38] delineated the first fuzzy clustering 
method "based on objective functions, and foreshadowed the 
usefulness of information measures (entropy) in the 
context of fuzzy sets. He further enlarged and illustrated 
his technique in references [39, 40, 41]. In reference [42] 
Dunn has defined the first fuzzy extension of the classical 
within group sum of squared errors (¥C00) objective func- 
tional. This algorithm is now known as ’Fuzzy ISO DATA’ 
process in literature after the hard ISODATA process of Ball 
and Hall [43] which is a special case of the former. 

Bezdek [44] has generalised this approach to yield an 
infinite family of fuzzy algorithms. Dunn [45, 46] has 
given further results on fuzzy ISODATA. In reference [47] 
Bezdek has given an interesting physical interpretation of 
fuzzy ISODATA process. Methods of clustering based on 
fuzzy graphs are still in their infancy, references [16,17, 
19,48,49,50] are seminal works in this direction. Some of 
the important applications of fuzzy clustering include 
plant classification [5l], medical diagnosis [52], and desigi 
retrieval, group layout and evolution of composite compo- 
nents [53] .. 
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2.4 Forecasting 

A large number of forecasting techniques have been 
developed during the last decades to fit the varied situations. 
¥heelwright[ 54 ] has suggested that these techniques can be 
grouped into two broad classes: quantitative techniques and 
qualitative or technological forecasting techniques. This 
classification generally reflects the extent to which a 
forecast can be based directly on historical data. Quantita- 
tive techniques are based on usage of past data whereas 
qualitative techniques are used when past data is not 
available. 

2.4.1 Quantitative Forecasting Methods 

In the area of quantitative forecasting methods we 
find a number of techniques whose common element is that the 
forecasts are based almost exclusively on past data. Some 
of the more widely used techniques in this class include 
moving averages, exponential smoothing, adaptive forecasting 
and regression analysis. Wheelwright and Mahridakis [54], 

Trux [55], Box and Jenkins [56] and Brown [57] have discussed 
these in detail. 

The simplest of all the forecasting methods are the 
moving average techniques. The exponential smoothing 
techniques are basically moving average techniques but they 
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represent all of the data and therefore are statistically 
stable. Exponential smoothing gives weigh tages to the data 
proportional to their age in an exponential fashion. Winter’s 
three parameter exponential smoothing model [58] and Brown’s 
single and multiple smoothing models [57] are significant 
contributions in. this area. These models cover regular 
demand with stable trend and seasonality. Other variations 
of these techniques take into account the irregularities in 
demand and the nature of the market. Porecasting of sporadic 
demand, strongly fluctuating demand, sporadic demand with 
homogeneous and heterogeneous markets and regular demand 
with heterogeneous markets are included in such variations. 
Trux [55] has described these cases in detail. 

The ability of an exponential smoothing model to track 
changes in the time series is dependent on the value of the 
smoothing constant. A relatively large value of the smooth- 
ing constant (close to unity) will cause the forecaster to 
respond quickly to changes in actual demand, reflecting a 
fraction of random changes in demand as well as actual shifts 
in the average demand. A small value of the constant 
(close to zero) will respond more slowly and smoothly. 
Therefore a fixed value of the smoothing constant in case of 
varying demand is likely to give erroneous demand estimates. 
Adaptive models overcome this objection to a great extent by 



continaally revising the value of the smoothing constant 
when changes in time series occur. Chow [59], Dudmsn [60], 
Trigg and Leach [6l] and Roberts and Reed [62] have contri- 
buted in this direction. Jain [63] has developed a hybrid 
forecaster using simple exponential smoothing and tracking 
signal model of Trigg and leach. This hybrid forecaster 
enjoys the strong points of both the forecasters as they 
are mutually compensating. The la.test addition to the 
family of adaptive forecasters is the adaptive filtering 
model of Wheelwright and Makridakis [54] which assigns 
weightages to past data -such that the mean squared forecast 
error is minimum. Rao and Shapiro [64] have developed a 
highly sophisticated adaptive forecaster using evolutionary 
spectra. 

Another class of sophisticated quantitative tech- 
niques is causal models. In these models the relevant 
causal relationships are expressed mathematically, leading 
indicators and regression models [ 57] and life cycle 
analysis [65] are the commonly used causal models. Some other 
complex quantitative models are Box and Jenkins method [56] 
econometric forecasting and input-output tables [54]. 

2.4.2 Qualitative Forecastins; Methods 

Qualitative forecasting techniques are used when no 
historical data is available. The aim of these methods is to 
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forecast changes in a basic pattern as well as the pattern 
itself. They are used primarily in two types of situations. 
First, to forecast when a given new process or product will 
become widely accepted and secondly for forecasting the new 
developments and discoveries that wjJ.l be made in a specific 
area. Because of the difficulty and cost of working wj.th 
qualitative methods of forecasting, they are generally 
applied only to long-term situations and to those of major 
importance to the fim. 

The mechanics of these techniques consists of inter- 
pretation of the related data by 'experts' who are supposed 
to ha.ve the best imderstanding of the situation and ability 
to predict changes of a similar nature. Market research [66], 
historical analogy [66], Logistics and S-curve approach [67] > 
time independent technological comparison [68], morphological 
search [69], Delphi technique [70, 71] and the relevance 
tree method [72] are some of the prominent qualitative 
forecasting methods. 

2.5 Inventory Control Techniques 

Inventory control has been one of the most favourite 
areas of interest for research workers and more analytical 
research has been directed towards inventory control than any 
other problem area in business and industry. Though the first 
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attempt dates back to 1912, significant progress has been 
made only after 1950 and v~JTy many realistic approaches to 
inventory control theory have been proposed. The inventory 
control problem has been tackled by analytical, graphical 
and simulation techniques and a good amount of research 
has been done in these fields. Recent application oriented 
research has been concerned over devising computationally 
efficient 'dynamic lot-sizing' algorithms for realistic 
situations in which the demand, rate of return and prices 
are variable. 

2.5*1 Analytical Techniques 

Newberry [73] has given a thorough account of the 
early developments in inventory control theory. In 1912 
Babcok for the first time tried to determine a mathematical 
basis for economic lot size by using a cubic equation. In 
1915, Harris developed an economic lot size formula which 
is almost identical to its present, foim. In 1918 Taft 
modified the economic lot size relationship to accoimt for 
the fact that manufacturing and sales period overlap each 
other since some of the production items are diverted to 
meet current orders without entering the inventory. T’ery 
few new developments were recorded between 1920 and 1931 with 
regard to the economic lot size formula. In 1931 Raymond 
suggested an approach to account for all the factors that 
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affect economic lot size and gave limits of possible lot 
sizes with little difference in unit costs. Around 1958, 

Sche rider, Erich and lessen tried to further generalise the 
problem by considering manufacturer's sales forecast, 
production capacity constraints, initial inventory, carrying 
charges and production costs. 

There were few, if any, major developments in the 
field of inventory control between 1938 and 1951. The 
attempts prior to 1951 considered deterministic problems 
but after 1951 researchers tried more realistic approaches by 
incorporating stochastic nature of demand in their models. 
Arrow et al. [74] in 1951 and Dvoretsky et al. [75,76,77] 
in 1952 and 1953 presented pioneering work in this direction. 
Bellman et al. [78] in 1955 analysed the dynamic inventory 
situation with linear ordering and penalty costs. In the 
subsequent developments the stochastic nature of lead time 
was also incorporated which made it necessary to consider 
the determination of reorder point together with order 
quantity. Since the reorder level requires stock to be 
reviewed to triggei a new replenishment order, and this can 
be done either continuously or at fixed interval of time, 
the further developments branched off based on the review 
procedure. ■ 


23 


The essence of continuous review models is to monitor 
the stock continuously as the demand occurs. This involves 
the determination of a reorder level, s, and the order 
quantity, Q, giving rise to two systems? the (s, S) and 
(s, Q) systems. In the (s, S) system an order, Q = S - s, 
is placed when the stock depletes to reorder level. Inventors?", 
therefore, fluctuates between two levels 'S' and 's'. In the 
(s, Q) system, a fixed quantity 'Q' is ordered when the 
inventory falls below the lower level 's*. 

The periodic review models are similar to continuous 
review models except that instead of monitoring the inventory 
status continuously, the review is carried out only after a 
time period 'T'. When 'T' is very small, the periodic 
review system tends to become identical to the continuous 
review system. 

A hybrid system, which is a combination of the above 
mentioned systems, also exists. This system is known as 
(s, S, T) or (s, Q, T) system. In this system, the replenish- 
ment is provided at every period 'T', but if in between the 
review periods the stock declines to 's', an order is placed 
to replenish the stock. 

ftirther literature is presented in the light of the 


above stated, policies. 
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Periodic Review Models 

The fixed cycle policy has received very little attention 
from the researchers because, theoretically, if the review- 
costs are ignored, all other things being equal, a continuous 
revie-t/f system always has a lower average annual cost than a 
periodic review system. However, this system sti3-l exists 
because it is not always feasible to make continuous re-views 
for low cost high volume products. 

Vazsonyi [79] has given an excellent analytical 
discussion for determining the maximum level, S, when the 
length of review period is given and demand is variable. 

Whitin [80] determined the optimal length of review period 
when the dem^d is variable. He calculated the upper level, 

S, by adding the expected demand during the length of the 
period to the safety stock required for any desired level of 
protection. Naddler [8l] and Churchman [82] have presented 
some simple models with variable demand which determine 
optimal value for the maximum level, S. Iglehart [85] 
investigated both back-order and lost sales cases and 
incorporated a variety of demand distributions. Yaspan [84] 
developed a model where the demand is a function of both the 
time required to fill the order and the nominal value. This 
is a realistic assumption because customer is less likely to 
order in future if he has to wait for back orders to be 


satisfied . 
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Continuous Review Models 

Continuous review models have received a very wide 
attention. Arrow, Harris and Marschak are the pioneers in 
this field. In their classical work [74], they considered 
fixed models, uncertain models, static models, dynamic 
models and also derived a procedure for obtaining optimal 
values of ’s’ and 'S’ in the (s, S) policy when the length 
of period is specified. Dvoretzky, Kiefer and Wolfowitz [75, 
76,77] have proposed three highly theoretical inventory 
control models. The first model is concemed with determining 
the optimal ordering policy when tho demand distribution is 
known? the second model deals with the optimal ordering 
policy when demand distribution is unknown, and the third 
model establishes the necessary and sufficient conditions 
for the validity of the (s, S) policy. *111 these models are 
highly mathematical and rather abstract and are not suitable 
for direct application purposes. 

In the literature, a few models which account for the 
safety stock required for varying levels of protection have 
also been reported, Yazsonyi [79] deteimined the safety 
stock for normally distributed demand, while Whi tin [80] 
tackled the problem for specified service level and 
Poisson demand. Buchan [2] calculated safety stock for 
specified service levels and for unifom, noimal and 
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exponential demands. Llewellyn [85] gave a cost model to 
calculate the lot size to he produced in order to compensate 
for scrap losses in manufacturing. 

The foregoing works give a general approach to 
inventory problems. In recent ye-Ts, different models 
for specific problems have been developed. lampkin [36] 
solved the problem for spare parts required for aircraft 
overhauling with Poisson demand and lead time with gamma 
distribution. A good number of models have been devoted to 
blood bank inventory control [87]. After the publication of 
perhaps the fiaist paper by Klein and Eosenberg [88] on peri- 
shable inventories in I960, it seems to be one of the live 
areas of current inventory research. Two recent papers in 
this direction are by Pries [89] and Nahmias [90] . 

In somewhat unconventional approaches, the inventory 
problems have been solved by queuing theory. Cohen [91] 
treated demand as an a,rrival process and replenishment as 
a service facility. On the contrary, Kaiush [92] treated 
demand as a service facility and replenishment as an arrival 
process. Decision criteria were server busy time, customer 
balking rate, waiting time for replenishment, etc. 

At present, the research work on inventory problems 
is being carried out at many different levels. At one 



27 


extreme, a considerable amount of work is concerned 
strictly with practical applications, while £it the other 
extreme, work is being done on the abstract mathematical 
properties of inventory models without any regard for 
possible practical applications. The practical situations 
demand simple and computationally efficient methods end the 
graphical and dynamic lot sizing techniques are oriented 
towards this. 

2.5.2 Graphic al Techniques 

Often mathematically sophisticated models prove to be 
futile in the face of the crudeness of available data. This 
difficulty has motivated researchers to carry out applica- 
tion oriented work including graphical, semi -graphical and 
empirical methods. Among the notable works are those by 
Herron [95], Eckey et al. [94], Hadley and Whitin [95] and 
Buckland [96]. 

Herron, using graphical and algebric techniques, 
has attempted the problem with stochastic demand and lead 
time by considering the cases of stockout penalty propor- 
tional to the number of stockout occasions and to the 
annual number of units out of stock. 

Sekey et al. give a semigraphical method for Poisson 
demand and lead time. Hadley and Whitin suggested an 
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iterative scheme for lost sales with stochastic lead times 
■and demand. 3uckla,nd developed a nomograph for normal 
demand and lead time and specified service level. 

2.5.3 dynamic Models 

Another class of practice oriented inventory models 
is dynamic lot sizing procedures 4 fhis type of procedure 
is called dynamic because it allows the time dependent 
pattern of orders and the rate of return on capital, as 
well as the pa.ttern of the unit price dependent upon the 
quantity, to be tehen directly into calculation as variables, 
whereas in the classical lot-size formula they are all cons- 
tant. The basic concepts of the method are as follows. 

In order to find the lot size for minimum total cost, 
we start from an order quantity whose total cost has been 
worked out. This order quant it j/- is then systematically 
adjusted and changes in cost arising from the adjustments 
are assessed. We then proceed in the direction of falling 
cost until the desired minimum cost is found. 

The optimum lot size is therefore discovered by 
means of systematic tests in contrast to the classical 
lot-size calcxfLation, which establishes the minimum cost 
by analytical methods which art- foasible only in limited 


c i r e urns t anc e s . 
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The Part Period Algorithm (PPA) developed hy 
Dg Mrtteis [5] in 1968, is the first effort in this direction. 
Its performance compares well with that of the Wagner - 
Whitin algorithm [97] which is based on dynamic programming. 
Though the Wagner-Whitin algorithm givc-s the best results and is 
used as a standard to compare the performance of various 
dynamic lot sizing procedures, it entails high computational 
burden and is too complicated to be understood by the average 
system user. The Part Period Algorithm overcomes these 
drawbacks to a great extent and the inventory cost deter- 
mined by it is only very slightly more than that from the 
Wagner-Whitin ■ algorithm. 

However, PPA is also not free from objections. 

Orlicky [98], who has presented an excellent review of 
various lot sizing techniques, has indicated some basic 
flaws in PPA. The PPA assumes that inventory depletion 
occurs at the beginning of each period, and a portion of 
each order equal to the first period requirement is consumed 
iminediately upon arinval in stock and thus incurs no inventoiy 
carrying charge. The inventoiy carrying charge is thus 
computed on the basis of this assumption rather than on 
average inventories in each period. Orlicky has further 
criticised the ‘look ahead‘ and 'look back' correction routines 
of PPA and has shown with examples that these special features 
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of PPA do not improve its effectiveness but on the 
conterary are dubious. 

Another work in this direction is by Pierce [99]. 

He considered a multi-item economic lot-sizing problem in 
which the ordering policies for individual items are 
interdependent as a consequence of joint ordering or pro- 
duction set up cost. The problem is to determine an 
optimal ordering plan in which the sum of the costs of 
carrying inventories and the costs of ordering are minimised 
and in which the known demands for each item in -each time 
period are satisfied. He has presented two algorithms 
to solve this problem. The first one is a direct algorithm 
which yields periodic solutions and applies to problems in 
which demand occurs uniformly over time. The second algo- 
rithm is based on dynamic programming and is an extension of 
Wagner-Whitin algorithm. 

Truz [55] has also discussed several dynamic lot 
sizing techniques which are very usefuil from the view point 
of computerised inventory control systems design. 

2.5.4 Simulation Models 

The analytical methods, generally, are v^^ry compli- 
cated if the demand and lead time both are stochastic in 
nature. The simulation of inventory situations on digital 
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computers has become popular due to great speed of modern 
computers and sometimes it is the only suitable method to 
solve complicated inventory problems. 

Quite a few developments have been reported in the 
literature regarding simulation of inventory situations. 
McMillan [100] gave a Monte-Carlo simulation study of a 
stochastic inventor;'?' problem. Cohen [91] formulated a 
multilocation model for a central store with branch ware- 
houses. Gross and Soriano [lOl] carried out a simulation 
study comparing the res'ulting inventory reduction obtained 
■when sealift is replaced by airlift. Lewis [102] simulated 
on a small hybrid computer industrial inventory situations 
in which demand a^nd lead time are both stochastic. 

2.6 Computerised Inventory Systems 

In previous sections a number of inventory control 
techniques were discussed "which reduce the managerial 
decisions at item level. This, however, expands the informa- 
tion for a single item to such an extent that it becomes 
almost impossible to carryout computations manually with 
speed and store the necessaiy variables and parameters for 
thousands of items in a perpetual inventory system, ilso, 
time is seldom available to forecast demand and to review the 
stock status for each item. Therefore, the help of computer 
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is sought to hsjidle the enormous amount of data at high speed. 
In past many a.ttompts have been mace to computerise a part 
or whole of the inventory system. Some of them are presented 
here in brief. 

Blondin [103] computerised an A B C system of inventory 
control in 1964 using IBM 7090 and IBM 357 data transmission 
equipment. The A and E items were tightly controlled and 
the key problem was receiving upto the minute records of 
stock balances, particiolarly for A and B items. 

A similar work has been published by G-eschelin [104] 
regarding computerised inventory control used for Chrysler 
car assembly group comprising of seven assembly plants. 

A third generation IBM computer was used. The method of 
operation wtcs based on Chrysler' s concept of management by 
'exceptions' such as 'shortages' and ' ove rshipmonts ' , The 
exceptions were considered to be the cause of delay and, 
therefore, required actions. 

Bext work in computerised stock keeping is an IBM 
Hospital Inventory Program [105] where a perpetual, inventory 
system provides for a continuous record of supplies, 
receipts and stock on hand. Items requiring majiagement 
attention are listed on exception reports. On a periodic 
basis transactions posted to the inventoiy record are 
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transferred to general ledger account which reflects the 
aggregate of receipts and issues posted to individual 
inventory records- 

In the foregoing wrks the computer effort was 
limited only to stock keeping and obtaining some useful 
reports. Further attempts have been made in the direction 
of greater automation to cover various management policies 
based on scientific inventory control techniques. One such 
effort has been in handling problems by LIMIT-Technique [106 J. 
The LIMIT program makes it possible to obtain the most 
economical lot sizes for all items in a selected group of 
items whore there is a limitation on the number of orders 
to be handled. 

Reed and Stanley [107] have designed a procedure for 
controlling general hospital inventories. They have used 
iBC invent or3/ control in which demands are forecasted 
for ii and B items by simple exponential smoothing. They have, 
however, not tried to find lot sizes optimally but the 
orders equal to the next period demand forecast are placed 
in each p^^riod. They havw given some indications on 
computerising’ the procedure. 
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Buchan [2] in his case study named 'Warehouse 
Inventory Management for a Drug Wholesaler*, has given the 
objectives of the system design and outlined the procedure 
based on Selective Inventory Management Concepts (SIM). 

The SIM includes the consistency in the ordering of inven- 
tory items by the routine use of ma.thomatical inventory 
formulae and selective allocation of available control 
effort to ea.ch item according to its importance. 

In 1969, Barret [108] discussed automatic inventory 
control techniques which include forecasting system, calcu- 
lation of order quantities for given discounts, allocating 
order quantities in a group order, setting, safety stocks, 
etc. He also discussed the advantages and disadvantages of 
various control systems in the light of computer application. 

The effort of Trux is of immense practical importance 
in the context of computerised purchasing and stock control. 
In his seminar papers, which were later published as a 
book [55], presented at the G-erman Institute of Industrial 
Management, Frankfurt, he elucidated various algorithms, 
procedural organisations and programs details regarding data 
processing for purchasing and stock control and materials 
management information system. 
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Next important work in th.is area is of Fetter [109] 
who has given an integrated computer oriented inventory 
management system. This work presents a centrally managed 
multi-item inventory system using either continuous or 
periodic review systems. He has listed out the initial inputs 
and computer opera.tions which include stock- reviews and 
calculation of reorder level and order quantity. 

In a recent book Orlicky [98] has discixssed the design 
and implementation al aspects of computerised material require- 
ments planning systems. Eecently Piiri [llO] has developed a 
computerised ordering system for drug distribution business. 

2. 7 Summary 

The review of forecasting and inventory control theory 
has shown that these problems can be solved by using a variety 
of approaches. The forecasting problem alone can be tackled 
by a number of techniques with va.rying degree of sophisti- 
cation. The systems designer is further stranded by the vast 
number of models available in inventor;/ theory. However, 
despite the wide spread use of selective inventory management 
concept, no effort has so far been made at evolving a general 
strategy which ca.n classify items in terms of their relative 
importance taking all relevant factors into consideration. 
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Tliough a vast amount of research has been carried 
out in inventorjr control theory, published papers rarely 
devote enough attention to the need for a linh between 
forocc^suing and reorder policy. further, we find very 
few publications devoted to hospital medicine inventory 
control problem which is quite different from the generah 
inventory control problems. 

In view of the above, we in this work have tried 
to introduce a general approach to selective inventoiy 
anagement using fuzzy sets and an integrated forecasting and 
and inventory control system has been designed to meet the 
requirements of medicine inventory control situation in 
hospitals. In the next chapter the details regarding the 
approach used in this work will be described. 



CHAPTER III 


PROBLEM STATEMENT AND METHODOLOGY 

The introduction and literature revie-w have already 
given an overview of the problem. In this chapter the 
problem and the solution methodology will be discussed in 
detail. The chapter comprises of three sections dealing 
with the item classification scheme, forecasting procedure 
and ordering ioolicy for medicine inventory. It is perhaps 
for the first time that a classification scheme based on 
the recently developed theory of ’Euzzy Sets' is being 
used in Inventory Analysis. Being a new approach, a major 
portion of this chapter has been devoted to describe the 
classification scheme. It is felt that the forecasting 
and lot-sizing procedures are relatively familiar and there- 
fore they are not dealt with in great detail. 

5.1 S t at ement of the Pro blem 

A general statement of the problem will bej ’Design 
a practical, efficient and readil3?' adaptable procedure for 
mana.ging medicine inventory in a medium-sized community 
hospital’ . 
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The first logical step towards designing any system 
is to specify its objectives and requirements. Before this 
is taken up, it would be worthwhile to consider what makes a 
system effective and efficient. Is the success attributable 
to the skill exercised in design, or to the caliber of 
management during ihe operation stage? Successful systems 
are often characterised by their simplicity, flexibility, 
reliability, economy, and acceptability. V7hile simplicity, 
flexibility and reliability tend to be function of design, 
economy and acceptability pertain to both design and operation. 
Though it may not be always possible to achieve this, the 
designer should aim at striking an optimal balance between 
these factors to meet the inventory management objectives. 

3*2 I nventory Managemen t Objectives 

The objective of any inventory management system is 
to provide the required goods, in the ri^t quantity, at the 
right place and at the right time, at minimum overall cost, 
further consideration will show that this essentially boils 
down to answering the following questions J 

(1) ¥hat should be ordered? 

(2) When should it be ordered? 

(5) How much of it should be orde3?ed? 
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All invent ory situations are alike so fa.r as the 
shove questions are concerned. However, each inventory 
problem will differ somewhat in the specific use of quanti- 
tative methods and the level of sophistication required 
in employing decision models. The supporting facts snd 
forecasts will depend upon the unique characteristics of 
the situation under consideration. 

The solutions are likely to be situation and problem 
oriented and this work is no exception. Though the bs,sic 
ideas and methodology described here can be applied to many 
inventoiy situations, what follows has a bias towards the 
problems of the situation studied in this work. Having 
seen what is expected from an inventory management system 
in general, we may now proceed to examine the hospital medicine 
inventory situation. 

3.3 The Hospital Medicine Inventory Control Problem 

To be successful, a management control system has to 
be compatible with the specific problem situation. The 
important characteristics of the hospital, medicine inventory 
anagement problem may be summarised as follows: 

i) There is a wide variety of medicines. 

ii) jPor some ailments, if the required medicine is not 
available, substitution by suitable medicines from 
stock is permissible. 
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(iii) In general, the requirements for medicines are seasonal 
with a definite trend. 

(iv) Some of the medicines are frequently in short supply 
in maxket. 

(v) The medicines ha.ve varying degrees of consequential 
importance in terms of probable loss due to non- 
availability in emergency. 

(vi) A small percentage of medicines tie up most of the ' 
capital. 

(vii) High service levels are recommended by hospital 
administra-tion as nonavailability of & medicine may 
be the question of life and death for a patient. 

(viii) The procurement lead time is variable for most 
medicines. 

(ix) Due to government control, the prices can be consi- 
dered stable for a reasonable planning period. 

(x) Special storing methods are required for some of the 
medicines. 

(xi) A good number of medicines, particularly injections, 
expire in due course of time. 

(xii) The present hospital regulations make it imperative 
to maintain stock recox-ds for all medicines. 

Though the above list may not be exhaustive, it 

nevertheless, gives a fair idea of the situation. The first 
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problem, as we have already seen, is how to deal with an 
inventory control situation when thousands of medicines are 
used. ¥e have also to decide the level of mathematical 
sophistication to be employed in the inventory control 
system so that an optimal comproiaise is sought between the 
system operating cost and the savings due to it. 

3 . 4 Selective Inventory Management Gonce-pt 

Selective Inventoary Management (SIM) deals with segre- 
gating inventory items into areas of relative importance. It, 
thereby, highlights the items for which the potential savings 
appear to be greatest and which should be controlled with 
utmost care. Only a careful and complete selectivity analysis 
in advance of any attempt at quantitative economic analysis will 
determine the extent to which the mathematical methods of 
inventory control should be applied. In such an analysis, 
each item in inventory should be ranked for importance in 
terms of pertinent characteristics. The procedure will vary 
from one inventory situation to another, it being necessary 
in some instances to consider all items to be of equal impor- 
tance. An example would be a production or assembly line 
where service level is important to keep the line operating. 
Where selectivity is warranted, which is the usual case, 
classification of items on the basis of relative importance 
can be ecnfirmed by an evaluation using an appropriate 
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mathematical model. The most commonly practiced, or perhaps 
the only known method for selective inventoiy management is 
the annual rupee usage based ABO analysis. 

3.4.1 Traditional ABC Analysis 

The basic idea behind the rupee usage based ABC 
item classification scheme is to apply the bulk of the 
planning and control resources to A-itcms, ’where the 
money is’ , at the expense of other classes which have demons 
trably much less effect on the overall inventory investment. 
In a typical ABC classification the top 20 percent, middle 
50 percent and the last 50 percent of the items may account 
for 80 percent, 15 percent and 5 percent of the total 
inventory investment respectively. However, there is 
nothing sacred about these figures and the proportions may 
significantly vaiy from situation to situation [7]. The 
steps required to perform an ABC analysis ares 

(1) Calculate the annual usage in units for each item. 

(2) hxtend the annual usage in imits by multiplying 

it with unit cost to get the annual usage in rupees 
for each item. 

(3) Bank the items from highest to lowest annual, rupee 
usage and assign categories. 

This procedure has been computerised in the present 
work and the results have been compared with those obtained 
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from the proposed 'Fuzzy Inventory Classification Scheme 
(FIGS)'. 

3 • 4 . 2 The lleed for a How Classification Methodology 

It is unfortunate that during the past years this 
type of rupee usage based analysis has been indiscriminateljr 
applied with no regejrd for the needs of a particular situa- 
tion. A good inventory cla>,ssif ication system should recognise 
many other characteristics of the items in addition to their 
annual rupee usage. It would be absurd to conduct an ABC 
analysis and apply a rigid set of ordering rules based 
solely on this criterion. To cite an exo^mple, in a hospital 
medicine inventory control situation, judgemental decisions 
about the importance of an item only on this basis may lead 
to serious trouble. There are factors like life saving 
importance, substituta-bility, expiry time, storage require- 
ments, market availability, lead time, demand and lead time 
variability and so on, in addition to annual, rupee usage, 
which must bo given due weightage in evaluating the relative 
importance of an item. As we have already seen, the principle 
of selective inventory management does not say that the 
judgement be based only on annual rupee usage. It simply 
seeks to apply maximum control effort to the vital few and 
allocate the control and planning resources in accordance 
with the relative importance of items, Ttirther, the 
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assignment of a certain percentage of the items to any paj?ti- 
cular class is rather arbitrajry. Bvcn highly dissimilar 
items may he put together and subjected to the same control 
procedure eveh when they resemble each other only on one- 
count, that is annual rupee usage. This calls for the need 
to evolve a methodology which takes into account all the 
pertinent characteristics of the situation, together with 
the inherent imprecision involved in the term 'relative 
importance' , and scientifically classify items in various 
groups. To this end, a classification scheme based on the 
theory of 'T'uzzy Sets' has been developed during the course 
of this work. 

3. 5 The Tuzzy Sets Approach 

More often than not, the cla. 3 ses of objects encountered 
in the real physical world donot hav^ precisely defined 
criteria of membership. For example, the 'class of animals* 
clearly includes dogs, horses, cows etc. as its members and 
clearly excludes such objects as rocks, fluids, plants etc. 
However, such objects as starfish, bacteria, etc. have an 
ambiguous status with respect to the class of animals. The 
same kind of ambiguity arises in the case of the* class of 
beautiful woman ' or the 'class of tall men' . Though they 
do not constitute classes or sets J.n the usual mathematical 
sense of these terms, the fact remains that such imprecisely 
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defined ’classes’ play an important role in human thinking. 

Ihe concept of ’fuzzy sets' is fast getting recognised as a 
powerful tool to solve problems involving such fuzzy ideas. 

A similar situation is encountered in Selective 
Inventory Management where the items are to be classified 
on the basis of their relative importance. It is not easy 
to answer what should decide the relative importance of an 
item, what is meant by terms like 'most important', ’\mimportant’, 
and 'vital', and which items are ’important' and which are not. 
In some way the ’class of important items’ is not any less 
vague than the ’class of beautiful women’ or the ’class of tall 
men’. 

The notion of a ’fuzzy sot' provides a convenient point 
of departure for the construction of a conceptual frame work 
which parallels in many respects the frame work used in the 
case of ordinary sets, but is more general than the latter 
and may prove to have a much wider scope of applicability. 
Essentially, such a frame work provides a natural way of 
dealing with problems in which the source of imprecision 
is the absence of a sharply defined criterion of class member- 
ship. The theoiy of ’fuzzy sets’ is on the verge of entering 
into an era of rapid development and several algorithms have 
been developed which have the potential to resolve the problems 
of ’fuzzy classification*. 



One such, problem is ABC type inventory classification 
which requires tho formation of groups of items on the basis 
of their 'relative importance'. These types of problems come 
in the domain of cluster analysis which is a well known 
technique of grouping objects. The traditional cluster 
analysis produces disjoint clusters, that is, an object gets 
allotted to one and only one cluster. Such strict classifi- 
cation is inappropriate, in most real life problems, like the 
one wo are faced with, because, of the imprecision or fuzziness 
present in the system. Fuzzy clusters account for this by 
assigning to every object in the cluster a me.Tnbership grade 
between zero and one. The classification system developed in 
this. work is based on the concept of fuzzy clusters and it 
successfully tackles the problem of forming . item groups on 
the basis of relative importance. 

Ve shall review now some important concepts in the 
theory of fuzzy sets and then describe how they have been 
eventually applied to solve our problem. 

5 . 5 • 1 Fuzzy Sets and Fuzzy Clusters 

The objective of cluster analysis is to group objects 
into clusters such that the elements within a cluster have a 
high degree of natural association among themselves while the 
clusters are relatively distinct from one another. With 
full faith in Pareto’s law of ’vital fevr’ , we resort to 
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cluster analysis to study and reveal the various importance 
groups present in. hospital medicine inventoiy. Our approach 
will he to first review some definitions, traditional cluster 
analjrsis, fuzzy clustering and its advantages, and then 
describe the fuzzy inventory classification scheme (FIC3) 
developed in this work. The review is based on the paper 
by Bezdek [ 22] . 


let, XCZR^, be a space £sf points (objects) with a 

generic element of X denoted by x. Thus, X =<Xps' •••> — n^' 

V/e shall call each x- = (x.,, x. ) a ’feature vector’ 

—1 il is 

whose j-th component, x. .= , represents j-th feature of the 
individual represented by Xj_. 


Any 'fuzzy set (class)’ A in X is characterised by a 
’membership (characteristic) function u.(x)’ which assigns to 
each object in X a real number in the interval [0, l]. The 
value of at x represents the ’grade of membership’ 

of X in A. 


Conventional cluster analysis, also known as 'hard 
cluster analysis’, sdeals with partitioning X into (c) non- 
empty subsets (clusters). The classification strategy is 
such that the partitioning is optimal with respect to some 
well defined classification criterion. Usually 2 < c ^ nj 
where (n) is the number of feature vectors in X. later we 
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shall see how the ’importance' of an item in inventory can 
be modelled into a feature vector and how the set of these 
'importance vectors' can be partitioned through clustering 
algorithms in order to find various distinct but natural 
importance groups present in it. 

The basic problem with all conventional clustering 
techniques is that even when wojdiing within the same metrical 
framework on X, they may yield different partitionings of X[4 2]. 
The inherent noniiniqueness of conventional clustering suggests 
opting for a fuzzy sets approach since fuzzy membership 
values incorporate a reflection of non-statistic al uncertainty 
about the structure in the data [22]. Further, the fuzzy 
characteristic functions portray the mathematical inter- 
relationships of the individuals in X, better than the hard 
partitions. 

Let be the vector space of all real ( cxn) metrices 

with the following inner product, induced norm and norm metric 
defined in it. 

^A, = trace (AB^) (3.1a) 

(a, a) = If A fl^ (3.1b) 

d (A, B) = }{ A - 3/1 (3.1c) I 

'V' A, B e the superscript (t) denotes the transpose 

operation. If, any U e satisfies the following conditions, 

we call it a 'hard c-partitioning of X'. Let u^^^ denote the 
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ki-th. element of U. 

^ki l<i£.n (3.2a) 

^ u, - = 1 j 'V' l<i<n (3.21) 

k=l " 

n 

21 > 0 ? ^ 1 ^ ^ C (3.2c) 

The elements of the j-th row of Uj = (u^ 2 _>***» 

the values of any characteristic function on X, say 
UjJ X—>{o, 1^ . The conditions, (3.2a) and (3.2h), together 
require disjoint intersections and collective exhaustion of 
X. Condition (3.2c) requires that each of the c -partitioning 
subsets is non-empty. The 'hard c -partitioning space, , 
associated with which is also the basic model underlying 
most hard clustering methods, can be defined as: 

= -fueY^ : U satisfies conditions ( 3. 2) f (3« 
c L cn J 

As we have already pointed out the conventional cluster 
ing methods are handicapped by their inability to partition 
X uniquely because the resulting partitions are dependent 
upon the structure of the data as well as the clustering 
algorithm. Hence we shall resort to fuzzy clustering to 
overcome this predicament. 
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Moreover, in real life we are frequently faced with 
the problem of classifying objects which may have mixed 
attributes. It might be inappropriate to assign an object 
to a particular group as it may have a certain degree of 
belongingness to every group. Conventional clustering 
models fail to accurately represent situations involving 
such *f\izzy uncertainty* in the structure of the data. 
However, the fuzzy model provides for such situations by 
allowing the individuals to have grades of membership in 
various groups. This can be achieved by simply extending 
the range of characteristic functions on X to the closed 
interval [0, l]. Thus u^^ : X-->[0, l] is a ’fuzzy subset 
(cluster)’ in X and the values Uj^j^ = called the 

’grades of membership* of the objects x^^' s in a fuzzy 
cluster Ujj.. We call, U e a ’fnzzy c~partition of X', 

provided it satisfies the following conditions: 

Uj^i e [0 , l] j V 1 < k < c ; 1 < i < n {3.4a) 

c 

u, - = 1 I 'V' 1 £ i < n (3.4b) 

k=l ■ 

J: u, - > 0 J ^ 1 < k < c (3.4c) 

i=l ^ 

The j-th row of U, denoted by = ^^jl’ * ’ * * ^jn^ ’ gives 

the membership values corresponding to the fuzzy cluster 
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! 2-^[0, l]. The ’space of the fuzzy c-partitions; 
associated with X' can he defined asJ 

e : U satisfies conditions (5.4) J (5*5) 


It has been shown that is, in fact, the convez 

hull of M in V „ (Bezdek f44l) and is also compact. These 
properties have helped in the mathematical analysis of 
clustering schemes defined on [22]. 

By incorporating fuzziness, the individuals in X can 
now have full membership in a sin^e cluster, or can claim 
partial membership in two or more fuzzy clusters. The condi- 
tion (5.4a) allows the various clusters to have overlapping 
intersections. Condition (5.4b) requires each Xj_ "to have 
full membership in X and (5.4c) requires that none of the 
fuzzy clusters is empty. The provision for an item to claim 
membership in several clusters simultaneously, is particu- 
larly helpful for the classification problems like the one 
we are dealing with. It enables us to figure out the natural 
pattern of ’importance distribution’ in the feature vector 
data of the medicines and also tells us how strongHy a panti 
cular medicine belongs to a particular ’importance group’ . 
This latter information will be helpful in finally ass idling 
an item to a particular importance group. 


i. . . ■ • 

C£.N s < - 


i'UK 

I:- RA^Y 

508«fr 
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Now the question hefoiTe us is, given liow can the 

fuzzy c-partitions of X he found? One nice way to optimally 
partition X is via the family of generalised WGSS (within 
group sum of squared) error objective functionals defined in 
[44]. On the cartesian product of with we define 
J) m e [l,oC) as: 


Jm ''f- !'> 


n 

i=l 


Jr (n- 
k=l 


•ki 




(5,6) 


where U e V (3^, . . . e and = (^^,..., 3^^) e R® 

for 1^ k ^ c, and jj . | is any norm on ’feature space’. The 
(c) vectors comprising V are presumed to have features 
prototypical of vectors in X and have a high affinity for 
membership in the respective fuzzy clusters • They are 

called ’centroids’ or 'cluster centers’ of their respective 
fuzzy clusters. The norm || * jj ’ which is a measure of 
similarity in (3.6) compares members of the data to each 
other indirectly via distance between them and the cluster 
centers. 

Optimal fuzzy c-partitions of X are defined as 

✓s A. 

part of pair (U, Y) which solves the problem. 

Minimize ; 

X) over 0 R®® 


(3.7) 
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which, is also called the gene re.lised minimum variance parti- 
tioning problem associated with X. It can be seen that for 
^ Jm = Jl-Y- 1 < Cc»C » reduces to the classical 

WG-SS criterion function. Detailed discussions on the merits 


and drawbacks of and its extension can be found in 
Wishart [111], Ling [112], Dunn [42, 45], and Bezdek [21,51]. 
Partitions obtained as part of solutions for (3.7) have been 
shown to be related to a well defined type of hard, compact, 
well separated (0¥S) clusters for X in [42] . The necessary 
conditions for solutions of (3.7) are derived in [44] for the 
class of fxnactionals in (3.6) whose norms are differentiable 
(e.g. , inner product induced norms). It has been shown that 


for m )> 1 and Xj_ 
order for (U, Y) 


to be a local 


, the necessaiy conditions 
solution for (3.7) are: 


in 




l<-k^c j 


(3.8a) 
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s (0, 1) fY’l^.k^c; (3«8c) 

The complete details for m = 1 and the singular 
case for some 'i' and *k' may he found in [42,44]. 

for m = 1, the requirement (3.8a) is replaced hy a 
nearest neighhouihood assignment rule and U e M is 
necessarily hard, the cluster centers in (3.8h) are merely 
the centroids of the hard subsets in U, and the resultant 
algorithm is essentially the 'hard ISODATA' process of 
Ball and Hall [43]. m > 1, equations (3.8) define 

the conditions for the 'fuzzy ISODATA* algorithm which has 
been used in this work. The mayor drawback with hard 
ISODATA process is that it ma3?- show the presence of 
clusters even though there may not be any natural groups 
present in the data set. The fuizy IbODATA clustering 
process largely overcomes this difficulty as reported by 
Dunn [42, 44] . 

After this review of fuzzy sets and fuzzy clustering, 
we now proceed to describe how the nedicine inventory 
classification problem can be solved by abstraction 
of the medicines into feature vectors representing their 
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importance and liow this can be used to effect selective 
inventory management for improved economic control of 
inventory. 

5*5.2 Import an ce through Feature Vectors 

In order to reveal the natural classes present in the 
medicine inventory by Puzz3r ISODA'IA process, the relative 
importance of an item has to be expressed in a vector form. 
The first step in this direction will be to objectively spell 
out which features or charaeteri sties of an item contribute 
to its importance and how can they be measu3^d. A suggested 
list of features for the vectorial representation of 
importance of medicines is given below. As stated earlier, 
we shall use the notation x^^^. to represent the j-th charac- 
teristic of the i-th item. Por ease in comparison and 
clustering, the features have been defined in such a way 
that all values lie between 0 and 1. 

(1 ) Annual hupee usage Gomponent 

The first important characteristic is the annual 
rupee usage of an item. The annual rupee usage of an item 
is obtained by multiplying the annual consumption in units 
'by the unit cost of the item. An item with higher annual 
rupee usage is more important as it entails higher invest- 
ment. let, 
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= cost per unit of item ’i' , 

= average yearly consumption, of item ’ i* in units 5 
n = total numler of tyj^es of items in inventory, 

Xii= annual rupee usage component of item ’i* . 

Then , . 

C. A. 

X = - — . X = 1(1) n 

Max [ C i = 1 to n j 

( 2 ) Exp i ry Time Oomponent 

Very often problems arise when due to lack of attention 
not enough time is left before the medicines can be utilised 
in order to prevent the loss due to expiry. An item with 
shorter expiry time deserves more attention than an item 
with relatively larger expiry time. 

Let, 

= expiry time for item ' i* in months, 
x^2 = expiry time component of feature vector of 
item ’ i* , 

= ““ L®1 ! ^ i = lU) n j. 


max 


- 


^12 


j '"r i = 1(1) 


n 


s 


max 


Then, 
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( 3 ) L ife-Saving Importance Component 

This characteristic may he compared with 'shortage 
cost' in industrial and cormnercial estahlishments . Some 
medicines are such that if they are not available in time, 
the consequences can be serious to the extent of causing 
permanent impairment or even dea-th of the patient. On the 
other hand, there are several other medicines like tonics, 
vitamin tablets and the like, which. do not matter much if 
there is a stock-out as their nonavailability will not 
cause any serious damage. It is difficult to measure this 
characteristic and therefor^' the value^^j have to be assigned 
to this component only intuitively. 

Let this component be denoted by x^^ and we assign 

values: 

x^^ = 0 j if the item is unimportant, 

x^^ = 0.5 ; if some damage may possibly be caused 

due to its nonavailability, 

x^^ = 1.0 j if the item is very important. 

This subjective assignment of values is justifiable 
on the ground that the fuzzy ISOLATA algorithm used for 
clustering groups similar items together and if two items 
are of equal importance tlaeywill be put in the same group 
and will be controlled by the same procedure. 
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( 4 ) S ubstitutability Oom-ponent 

Most medicines cen be substituted by others in 
several situations. If substitutes are always readily 
available we may not like to give as much attention to an 
item as we will give to an item which has no substitute. 

This characteristic is also difficult to measure and the 
values have to be intuitively assigned. 

let represent the substitutability feature 

of item 'i*. ¥e assign the following values^ 

= 0 } if the item can always be substituted, 

= 0.5f if the item can be sometimes substituted, 
x^^ = 1.0*, if the item can never be substituted. 

A similar situation is present in departmental stores 
where for several items such as soaps, tooth-pastes, shaving 
creams etc, one brand may bo acceptable to the customer for 
another. 

( 5 ) Market i.vailability Component 

There may be some items which are difficult to 
procure from the open market when needed such as baby foods 
and some type of injections. This is particularly true for 
grocery stores where several items are always in short 
supply. This component can be measured using the availabi- 
lity ratio relationship given as. 
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ITumber of times item not available from supplier 


Total number of times purchase orders made 


"V" i = l(l) n 


where, 

= market availability component of item 'i'. 

If past records are difficult to search, a simple 

measure based on subjective judgement may be used. ¥e assign 

vailues to x-r- asj 
i5 

^i5 = 0 » if item is always available, 

^i5 ~ 0*5 f if the item is onl 3 / sometimes not 

available , 

^i5 ~ 1.0 5 if the item is difficiiLt to procure. 

# 

( 6 ) Stock-Out Eistoiv Component 

It is desirable to include an indicator of the 
efficiency of past and existing control procedures for deciding 
the importance of a medicine. If the earlier procedures 
were incompatible and have resulted in frequent stock-outs, 
the item needs to be promoted to the class of moare important 
items because it requires better control. The stock-out 
history component, the i-th item is expressed as, 

Humber of times the item found out of stock 

js: = — ^ , 

Total number of checks made 
^ i =1(1) n . 
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As an alternative the values 0, 0.5 and 1.0 may be assigned to 
this component for itemr which were practically never out of 
stock, sometimes out of stock and frequently out of stock in 
th sA o rd er . 

( ^ ) Storing; he t hod component 

Some of the medicines, particularly injections, require 
special storing methods like low temperature, protection from 
sun light and humidity, etc. The stocks of these medicines 
should not be allowed to pile up to the extent that it 
becomes difficult to accommodate them. 

Let this component be denoted by 'x^,^’. ¥e assign the 

values: 

x^,^ = 0 5 if no special care is needed in storing, 
x^rj = 0.5> if some special care has to bo exercised, 

X = 1 j if special care is required. 

(8 ) Storage Space Component 

Sometimes when there is a liraitation on the total 
space available for storing, particular care should be taken 
in maintaining stocks of items occupying too much space in 
order to accommodate other items. Th Ough this characteristic 
may not be of much consequence in many practical situations, 
yet it needs to be incorporated for the sake of generality. 
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let us denote this component by Then, 

Space occupied by item ’i' on the average 
Maximum of the spaces occupied by all items 

^ i = 1(1) n. 

One can determine x^g by determining the space requirements 
of the -various items, ihi alternative proposition is to group 
the items into three categories and assign values of 0, 0.5 
and 1.0 depending upon whether the storage space requirements 
are low, medium or high for the item. 

(9) Lead Time Component 


The procirrement lead time adds significantly to our 
feeling of importance for an item. Needless to say, items 
with higher lead time need greater attention, let, 


i9 


T. 


1 


T 

max 


Then, 



lead time component of item ’i’, 


lead time of item ’i’ , 
Max I Tj, ; ^ i =1(1) n^ . 


T. 

X 


T 

max 


j i = 1(1) n . 
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( 10 ) lead Time Variability Oompoiient 


If the procurement lead time for an item is highly 
varying it deserves more attention and larger safety stocks. 
Let, 


X. 


1 10 ~ lead time variability component of item ’i' , 


Orp = variance of lead time for item 


Then, 


^i,10 




Meo: j(a^j /r. l(l)n ^ 

K, ± 


; ^i = 1(1) n. 


If it is difficult to calculate 8.nd cr the values 0, 

0.5 and 1.0 may he assigned as an alternative to items with • 
low, medium or highly varying lead times. 

( 11 ) Demand Yariahility Component 


If the demand for an item is highly fluctuating, it 
needs a better forecasting method to account for all the 
factors due to which the variations occur and to track the 
changes. One of the measures may be the ratio 
of variance with respect to mean demand. Let, 

Xi 11 = demand variability component, 

D^ = mean demand of item ’i* per period, 
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a = variance in demand per period, 
i 


Then, 


X. -IT 

1,11 


E. / Dp 


Maz|.0^j. /D. jVi = 1(1) 

^ i J 


5 "V i = 1(1) n. 


(12) Quantity Discount Oomponent 


The bulk purchase quantity discount may vary from item 
to item and manufacturer to manufacturer depending upon their 
marketing policies. Items for which higher discounts are 
available deserve more attention and care should be taken 
to use these discounts to advantage. let, 


^i,12 


^i 


Then, 


quantity discount component for item *i’, 
percent quantity discount available for 
item ’ i' , 

cost per \mit of item 'i', 
annual consumption of item ’i’. 


X 


P.C.A. 

Ill 


i,12 


MaxfP.G.A. jY"! = 1(1) a 

1 X i 


1 


J^i = 1(1) n 


The above list of features may not be exhaustive but it gives 
an idea regarding how to go about defining and expressing the 
feeling of importance for a medicine through feature vector . 
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Total space occupied by 'DRISTlif' tablets is 'too much' 

Procurement lead time is 40 days with a variance of 
10 days and maximum lead time for any item is 100 days. 

Maximum value of variance to lead time ratio for any 
item is 0.5. 

Ratio of variance to mean demand =0.4, and the maximum 
value of this ratio for any item is 1.0. 

Maximum^ quantity discount available for 'ISISTilT' 
is 20 percent, and the maximum value of the factor 
^i^i'‘'u item is Rs. 3000/-. 

From the above data and the definition of various 

of 

feature vector components, the feature vector'l ' pRIbT/lT' iss 

^ = 1^0.25, 0.75, 0.5, 0, 0.5, 0,0. 5, 1.0,0. 4, 0.5, 0.4,0. 

After conversion of the information a,vailable about the 
medicines in stock into feature vectors, the next step is to 
classify the items into various importance groups. 

3.5.3 Olassification by Fuzzy ISORATA Process 

As we have already stated, conditions (3.8) together 
with convergence criterion and tie breaking lules form the 
basis of a process of simple iterative optimisation of 
defined in equation (3.6), called fuzzy ISODATA which has 
been used here to form clusters from the input of feature 


(x) 

(xi) 

(xii) 

( xiii) 

(xiv) 


vector data. 
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The algorithm starts with a set of *c' initial centroids, 
(i'he next section deals with the selection of these centroids). 
In every iteration the membership grades of each item with 
respect to every cluster are determined and the centroids 
are recalculated. The iterations are stopped when two 
successive sets of centroids are very close to each other. 

The steps in the algorithm are as follows: 


(1) Start with 'c' initial centroids y^ , j v . 

(2) Set i = 1. . 

(3) Consider the feature vector of item 'i*. 

(4) If is equal to any of the centroids, say y., compute 

a- J 

the membership grades of x^ with respect to the clusters 
by; 


u^(Xi) = 1 J for i = j , and 

u.(x. ) = 0 J for j 5 ^ i. 

J 

If is not equal to any of the centroids, then: 


(%) = 


c 

r 

1=1 


It^i - 

|!% - vpli 



J for V" k =1(1) c 


(3.9) 
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where Ij . | i 


i is a norm induced by the weighted inner produc 


’ ~) 


s 

on and 


2 2 

A = [ dlag ( -1~ -®- )] 

w 2 ^ w 


= diag ( “2 , 


P— M-1 
s 

Zi 

• • • 9 2 


), 


j cr _ variances of the feature vector 
s ^ 


components, 


w 


2 _f •.•>Wg normalised weightages given to 


various feature vector components, 
number of feature vector components. 


Upon substitution of the norm in (3.9), we get. 




X. .-V . 


2 . 

3=1 ^ 3 


(l/(m-l)) 


c 

H 

1=1 


s 


i/ r w . ( ) 

/ 3=i ^ 


j 


(— i — ) 


k = 1(1) c ^3 


(5) Set i = i + 1 

(6) If (i^ n) go to next step, otherwise go to step (3)j 


where 'n* is the number of feature vectors, 


lo) 
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(7) Calculate iiaw CGiL"troids by: 


n 


:?k = 


.£ [U4.(x. )]”* J. 


n 


(3.11) 


[“k {5i)] 

1=1 


m 


(8) If, 


,A 


. 11/2 


! -C ^ 5 forV" k = 1(1) c, stop 


D=1 


Otherwise put and go to step (2), where 

'e' is a small quantity whose value depends upon the 
accuracy desired. 


¥hen the process is complete, final centroids together 
with the membership grades of items in each cluster are 
printed out. As the algorithm treats all the elements in the 
feature vector as continuous variables, the centroids do not 
represent any real item. For a given set of weights, the 
fuzzy clustering process is repeated for several values of 
*c’ and 'm*. It should be noted that for any item the 
membership grades sum to unity over all the clusters. 

Figure (l) depicts the logic of Fuzzy ISODAIA process. 

At this moment two questions remain unanswered. Firatly how 
should the initial centroids be chosen and second what would 
be the best choice for Vc’ and 'm*. In the nezt section we 
will try to answer these questions. 
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5 • 5 • 4 Ohoosin^ the Initi al C e ntroids 

The fuzzy ISODATA is sensitive to the choice of 
initial centroids [53]. So far no absolutely satisfactory 
method has been reported for this purpose. Anderberg [30] 
has given several methods for choosing the initial centroids. 
.In this ■work, we have used a method suggested by Ball and 
Hall [7.13]. The m-oohca can be summarised in the following 
steps : 


( 1 ) 

( 2 ) 


(3) 

(4) 


Start with the set of feature vectors x-, , x . 

Fix the first centroid v^ as the mean feature vector 
whos'^' components are given by: 





n 


i ^ j 


1(1) s 


and set K = 1. 

Set i = 1. 

If is same as any of the previously fixed 'K’ 
centroids go to step (6); bthermse, calculate the 
distance ’d- ■ of the feature vector x^ from every 

centroid by: 



s 

= / 7 w . 

uti ^ 



} "V" P = 1 ( 1 ) D • 
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(5) If ^ © j p = 1(1) Kj , set Yjr+p ~ = K+1, 

and go to step (3). where 9 isaconveniently set quantity 
denoting that all initial centroids should be separa- 
ted from each other atleast by this distance. 

(6) If (i = n) and no new centroid could be found go to 
next step, otherwise set i = i+1 and go to Step (4)» 

(7) If (K = c) stop. If (E < c) chan^ © to (@ - >*) ; 

If (E > c) change 9 to (0 +,X) and go to Step (2), 
where X is a. small quantity. 

The centroids thus generated for various values of 'c' 
are taken as input to the fuzzy ISODATA, The next section deals 
with finding the optimal number of clusters (c) and the 
effect of exponent (m) on clustering. 

A logic flow diagram for the selection of initial 
centroids is given in Fig. (2). 

3 . 5.5 Measures of Partition Quality 

For the first time Bezdek [22j has suggested a method 
for systematically finding the number of clusters in X. 

Though so far only a limited number of experiments have been 
conducted with this method, we have found it quite useful 
during the course of our work. The method is based on the 
concept of ’ quality' associated with fuzzj’’ clusters. 








WMm 
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Irrespective of the procedure used for clustering 
we may have to answer q,uestions like: does the algorithmic 
output really correspond to structure in X? More generally, 
does it correlate with natural cluster:^ in the real model 
that X represents? In all but trivial instances [ 22j , the 
answer to these questions is easy, we simply do not knoxf. The 
only criterion capa.ble of evaluating algorithmic clusters 
is undisputable a priori labelling (i.e. pre-partitioned 
samples), in which case there is no problem to solve. 

The difficult’’’ arises because 'cluster’ becomes a 
meaningful term only in a well defined sense usually deriving 
from the na-tui^' of the data, and of course the model we 
choose to represent it with. It appears that the heuristic 
approach of finding an algorithm capable of reliable decom- 
position of large training samples from the process under 
study is the only feasible alternative when dealing with 
hard clustering techniques. On the other hand, the fuzzy 
model provides some added information about structure in 
X which seems useful in evaluating the relative merits of 
competing classifications independent of the clustering- 
procedure used to obtain them, Bezdek [22] has suggested 
two measures of evaluating the partition quality. 
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Tile first measure is via the ’partition coefficient’ 
Tq(U) which can he defined as follows: 

2 

J’c(U) = trace (Uu‘'^)/n = — 51iO — .^11 e (3.12) 

n ^ ^ 

The values of I'^(U) are inversely proportional to the average 
content of fuzzy intersections in U and it intuitively 
provides a measure of the uncertainty an algorithm 
experiences when assigning membership values to individuals in 
It seems plausible that an algorithm yielding fuzzy 
c -partitions with higher value of F„ is more successfully 
decomposing X than one whose partitions are inherently 
' fuzzier’ . 

A second approach to this question arises from 
similarities to the notion of uncertainty usually used 
in information theory. Probabilistic measures of uncertainty 
in connection with classification models have been inves- 
tigated by Bstabrook [ 116 ], Orloci [ 117 ], and Wallace and 
Boulton [ 118 ]. Buspini [ 38 ] and Deluca and Termini [ 13 ] 
have discussed the notion of entropy in the context of 
fuzzy sets. Though there are philosophical distinctions 
between probabilistic and fuzzy uncertainties , the formal 
similarity between them suggests exploiting sane mathematical 
properties known in information theory to describe statisti- 
cal uncertainty. The ’nozmalised cla3Sificatio^ entropy * 
in U, denoted by H^(U) , based on the Shannon’ s entropy 
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function, suggests using a functions 


H„(U) 


i 


n c 

r z 

i=l 1=1 


u. , 
xl 


lo 


7 


^^il^ j 


n 


(5.13) 


T^rith Uj_^ log (Ujj_^) =, 0, if u^^ = 0. 

Although the functional forms of and H are quite 

o c 

different, the relation hatween them can he shown to be [22] 
as follows. 


Fq (U) = 1 <=r^> = 0 <=r) U £ is hard (-3.14) 

^c c <=> U = [1/c] (3.15) 

Further, Bezdek [22] has proved that the upper and lower 
bounds on H^^(U) ares 

H^(U) •, 

1 - Fq(U) < ( 3^1^ )< 1 (c - (U)) (3.16) 

For a = e this reduces tos 

1 - F^ (U) < H^(U) < i (c - F^(U)) (3.17) 

Equations (3.14,15) suggest that the equi-membership 
partition IJ = [-] , i.e. u^ = ^ jY i, t * is the fuzziest or 
worst one can do? on the other hand, the ideal situation 
occurs when the structure in X is so distinct that a fuzzy 
algorithm recommends a hard c-partitioning of X. Maximising 
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different fixed points of fuzzy algorithms minimises 
the total content or overlap in pairwise fuzzy intersections? 
equivalently, minimising over the same choices maximises the 

information extracted from U. These measures serve as a relative 
indication of the uncertainty an algorithm experiences in 
trying to assign memherships to the vectors in X. 

Numerical evidence [ 22j indicates that 'H ’ is more 
sensitive them ’Nq' in ranking U's? this is due to the fact 
that the slope of the logarithmic cirrve on most of (0, l) is 
much steeper than that of a parabola. In general the cluster- 
ing strategy has been to try fuzzy ISO DATA for clustering 
for a given set of weights, over several values of 'C and 'm' 
and regard that c -partitioning of X as optimal which yields 
minimum value of ’H * . The scheme has been shown in Hg. (3). 

In the present investigation of medicine inventory both ; 
the measures of partition quality were tried. It was found 
that both of them Suggest the presence of three partitions. 

The detailed results regarding these are presented in 
Chapter IT. The first group consists of medicines which are 
important from all considerations, second group has medicines 
which have low annual rupee usage and high expiry time but 
are otherwise important and the third group comprises of 
unimportant medicines. For assigning the medicines to various 
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clusters, the 'maximum membership grade rule' has been followed. 
Some other rules which use a higher cut-off value are also 
available. However, to avoid misclassif ications it is better 
to present the membership grade distribution information as 
it is and then make suitable transfers within various classes 
if necessary. 

The next task before us is to select or devise 
suitable ordering schemes for various classes of medicines 
and the remaining part of the chapter is devoted to this. 

3.6 Selecting the Ordering System 

In inventory models we attempt to isolate variables 
and parameters and try to understand the behaviour of 
inventories, ■‘^'urther, we seek to find optimal values of 
these variables under a wide variety of conditions and 
different degrees of model complexity. Then the management 
selects an inventory control system for implementation based 
on a number of criteria, many of which conflict with each 
other. An important criterion in a system's selection is 
its adaptability to existing procedures and regulations. 

This is necessary to minimise the resistance encountered 
while changing over to the new system and the implement at ional 
costs and delays. Therefore a look into the existing proce- 
dures for controlling medicine inventory in hospitals is 

desirable for evolving the new system. 
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Tile present regalations in most hospitals ask for 
maintaining perpetual stock records for each item. 
Medicines are issued to compounders from stores, as and 
when needed, and the transactions are entered in a ledger. 
When the stock of an item has depleted to some pre-set 
reorder point, a replenishment order for a suitable quantity 
is placed. At present, in most of the hospitals, no 
scientific method is used to calculate the reorder level 
and reorder quantity and the decisions regarding these 
are taken mostly on the basis of intuition. 

Prom the above, it is clear that a continuous review 
reorder point system will be readily adaptable to present 
regulations. Here we have two alternatives, ¥e can either 
opt for a conventional fixed limit control system or choose 
an adaptive control limit system. This decision is not a 
difficult one to make as the superiority of adaptive limits 
for inventory control is now more than established. Hilon 
and ElmaJ-eh [114] on having simxilated about 200 cases over 
500 periods, leported 50 percent reduction in inventory 
Costs over fixed limit controls by using adaptive control 
proesdures. Trux [55] has also stated that 30—40 percent 
savings occur by using adaptive limits- 
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Thougii we have seen that the 'once and for all* 
determination of control parameters based on historical 
data is clearly unsatisfactory, the analytical treatment 
of iii8.theniatioal models incorporating the real life features 
of stochastic demand and lead time is too complicated to 
be useful for practical applications. Wagner et al. [I15j 
have demonstrated the computational difficulties involved 
in determining control limits in such cases. 

In view of the above, we have suggested a method here 
which is based on the works of Buchan and Koeningsberg [2] 
and BeMattis [3]. This method for computing adaptive control 
limits is simple and straight forward, and it takes into 
account seasonal fluctuations and trends in demand. It, 
however, does not seem to be economical for the class of 
unimportant* items and therefore fixed limit controls 
have been used for these items. Before the details of the 
procedures are described, it is desirable to have an overview 
of the system. The next section deals with this. 

3*7 An Overview of the S ystem 

The general procedure will be to first locate *c » 
initial eentriods by Ball and Hall's method. Starting with 
these initial centroids, we classify the items for several 
values Of *c* and 'm* and accept that classification which 
yial(js best values of the measures of quality, and . 
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liiis classification needs to le done only once in a year 
or when there are a large number of additions or deletions 
in the list of the items. • This information together with 
the demand, cost and lead time data is taken as input to 

the program for computing lot sizes and reorder points for 

specified serYice levels. The program uses different algori- 
thms for computing these depending upon the class to which 
the item belongs to and gives a detailed report for the 
purpose of control. The procedure has been shown in Fig, [4]. 

The reorder level should be noted on the same page of 

the ledger which is being used for keeping the stock record 

for the item. When the stock falls below the reorder level, 
then the program which works on independent ordering 
policies may be run for only those items which need to be 
reordered. 

Though in its present form the system is only semi- 
computerised, it can be extended for automatic data process- 
ing. The present card files can be replaced by storage 
devices on a computer and the transactions may be directly 
entered. The reorder level information can be stored side 
by side and whenever a transaction occurs, a computer check 
can be made to find whether an order needs to be placed. 

Having described the genersil logic of the system, 
we shall now elaborate different procedures used for 
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controlling the inventory of three classes of medicines. 

5*8 Qontrol_^Pro ^^ure. for » Class I» Medicines 

The first class consists of those medicines which are 
important from all considerations. Keeping ia view their 
high importance, an adaptive control limit procedure has 
been used. The procedure essentially consists of forecasting 
the monthly requirements for each item and at the same time 
calculating the total unit inventory cost. The calculations 
are terminated when the lot size for minimum unit total 
inventory cost has been determined. This is checked against 
the expiry time of the medicine. The forecast error is 
also accumulated side by side. This together with the lead 
time variance is used to find the variance of the demand 
during lead time. This information is then used to set the 
reorder point for specified service levels by using Gamp and 
Meidel^s extension of Tchebycheff ' s inequality. The control 
parameters are recalculated when a new order is placed for 
replenishinent by updating the forecast. 

3.8.1 The Forecasting Scheme 

The ordering decisions are based on the demand forecast 
for each item. The demand for medicines exhibits seasonality 
and an upward trend in general. The trend can be attributed 
to the population growth siitf-in^^easing' medical facilities. 
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Seasonality is due to th.o climatic changes. A forecasting 
method suggested by Buchan and Koenigsberg [2j matches with this 
situation and is economical to use for a large number of types 
of medicines. The following notation is used for describing the 
forecasting procedure. Bor an item =iS let, 







a 


^i 






= demand in period -y of year 'l^» 

= number of years for which past demand data 
is available, 

= number of periods In a 3^ear (H - 12 in this 
case), 

= demand forecast for period ’t', 

= smoothing constant, 

= trend per period in demand, 

= base series ratio for period • j , 

= smoothed ratio upto period '3', 

= base series value for period ’3’, 

= number of periods of current year for which 
demand data is available. 


The steps are: 

[Ij Compute the trend per period in the demand for the 
medicine by; 



evelop the seasonal base by eliminating the overall 
monthly trend from the average monthly demand 

series as follows; 


2 ^.. , J /I 


+ t(| - 0.5) - j ]T. I 


Gompute the smoothed ratios by: 


1(1) IT . 


D. 

" B. . ^ 'ik'. 

IJ 


~t 3 = 1(1) ITs 


1 f 


\,J-1 » ^ j = 2(1) B. 

^ic ^ smoothed ratio series is updated upto the 

^^j_c~th period of the current year. 

[4] The demand forecast pf the,lj^-th period ahead is given 


'i,(t+L) = h,t-l ®(t+L,-l) h ■'i 
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Tiio demand forecast is taken, as input to the lot 
sizing model which is descrihed next, 

3.8,2 Lot-Sizing Procedure 

A dynamic lot-sizing procedure has been used for 
the 'Class— I’ medicines to determine the order quantities. 

It is a modified version of the Part Period Algorithm [3]. 

The modif ica,tion is made in calculating the cost which is 
calculated here on the "basis of a,verage' inventory since no 
depletion is assumed immediately upon arrival of the medicine 
in stock. P'urther, in case of medicines the prices and 
carrying costs can be assumed constant over a considerable 
period of time and, therefore, the dynamic lot-sizing model 
used here assumes only the demand as variable. The following 
notation will be used in describing the models 


Por any medicine 'i', let. 


C. 

im 


\s 


Pic 

C. 

IS 

N 


total cost per unit if the lot covers 
requirements upto m'-th period, 
stock-holding cost per unit per period 
in Rs . , 

price per unit in Rs., 
carrying cost in percent, 
set-up cost per order in Rs. , 
number of periods per year. 
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F. . 

13 


= demand forecast for period ’ j’ in units, 

= expiry time for themodicine in periods, 

= annual demand of the medicine in units, 

= optimal lot-size in units. 

It can be noticed here that the stock-holding cost per unit 
per period is given by: 


h . 
1 

Ai 

Ql 


^i- / 

^is IOOTII” Rs. /period /unit. 

The lot-sizin,;', procedure can be summarised in the following 
steps ; 

(1) Oompute ~ forecast and find 


hi = 


' hs ths hi 

hi 


( 2) Set, m' =111' + 1 


( 3 ) Forecast the demand F^^^, for period ’m’*. 

(4) Compute the total cost per unit by: 


im' 


m' _ 

hs + h. hj(3-0.5) 

m' 
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(5) 


If, 



m’-l _ 

= Zl.. F. • and stop. 
j=l ■ 


= C. , ^ 

i,m’ -1 

If not so 


go 


to Step (6). 


in' -1 

(6) If, m' > take = 2II B . - and stop, 

j=l 

otherwise go to Step (2). 

After calculating the lot-size the next important step 
is to determine the reorder point for desired service level. 
This is described in next section. 

3.8.3 Oalculation of Reorder Point 


An elaborate pixDcedure has been used to determine 
the reorder level in case of the first group of medicines 
as they are important and tie up a considerable portion of 
the total investment in inventory. The procoduie essentially 
consists of calculating the demand during lead time and its 
variance by using the lead time variance and forecast 
error. By scientifically adding adequate buffer to the 
average demand during lead time, desired service levels 
can be expected to be achieved. 

The forecast error can be calculated at the same time 
when forecasts are being generated by using the following 
relationship: 
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cr 


2 

if 


- 51 . 

j=i 




ij 


- ■)' 



( 3 . 18 ) 


wh.e re , 



^ij 

= 


estimated demand of item 'i' for period 'j*? 
actual demand of item 'i* for period *j’, 

= forecast error for period *j's 
variance associated with the next period 


forecast, n’+l’ item 'i', 

= number of time periods used in variance 
determination. 


To be able to set reorder levels for various medicines for 
desired service levels, it is necessar3?' to know the lead 
time frequency distribution. It has been found that in 
hospitals, in general, the stock items may be purchased from 
several vendors and the method of purchase from the same 
vendor may be by either bid or direct purchase order. These , 
differences pose several problems in computing lead time 
parameters for the medicines. To overcome these restrictions 
total lead time has been broken up into three subelements and 
some assumptions have been made which are as follows; 
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Let the total lead time be denoted by and further 
let ^i2 ^±3 three eubelements , where, 

^il ~ time required to write and receive approval 

of the purchase requisition, 

time required to write the bid quotation 
sheet plus, the time taken by vendor to 
consider this bid, and the time for the 
hospital staff to wite the purchase order 
after accepting a bid, 

I'ij = the elapsed time from issuance of the purchase 
order until the units are ready on the shelf 
for internal distribution. 

¥c assume that 'Tj_2' same for all. 

medicines, i'heir values can be established from .past 
data of all items and periodic review by sampling can be 
made to prevent any error due to abrupt changes in proce- 
dure. Por nonbid items, time will be equal to 

purchase order writing time, sinco the sublement bid time 
will be zero. ¥e further assume that the three times are 
independent and, therefore the estimated mean value of the 
total lead time for an item can be taken as sum of the 
sampled mean values of %1* % 2 '^13* 

with a bar we have. 


Denoting averages 
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(3.19) 


and xrom. the theorem of linear coinhinations applied to 
sample values. 




+ 




( 3 . 20 ) 


Lxt the demand during the lead time of an item he denoted 

by the mean of the sample distribution by X^^ , and 

2 

the varia.nce by o . Then for the i-th itera 

1 


X. = 
1 




m' 


and 


( 3 . 21 ) 


= 4f ^i + ^ “ 

^i 11 1 m 


i ^2 . 2 

i 


(3.22) 


where is the lot size tha.t will meet the requirements of 


m' periods alieacl from now. 


Since both the demand distribution and lead time 
distribution are unimodel, it can be safely assumed that 
the product of the two is also unimodal. The reorder level 
can then be computed using Camp and Meidel’ s extension of 
Tchebycheff ' s inequality which states: 


P(X.> Z. + a ) < 


4.5 z| 


(3.23) 
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where , 


= probability of stockout 


= (l-service level ’p'. ) 

^ 10 


+ Z ■ O’ 
1 1 


= order point providing a proba- 
bility of 'p^q' of not having 
a stockout during replenishment; 

Pio ~ service level = probability of 

not being out of stock. 

The inequality can be rewritten as. 


Z. > 


V 4.5(l-Pio) 


and , 




X. + Z. o 
1 ^ 


(3.24) 


(3.25) 


where , 


ho 


= reorder point for medicine ’i' 


Thus for a prescribed service level, 'Pj^^’ > can be 
computed readily from relation (3.24) which upon substitution 
in equation (3.25) gives the reorder point This completes 

the determination of all control parameters for ’Glass-I’ 


medic ines . 
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3 • 9 Control Procedure for * Glass -I I' Medicines 

The second group is made up of those items which 
have a low annual rupee usage but are otherwise important. 
Due to this reason the same foreca,sting and lot sizing 
procedures have been used for this group as were used for 
the 'Class I' medicines. However, the procedure for 
calculating the reorder point for 'Class I' items is not 
economical for this class as the advantages gained from 
tighter control may be offset by the increased cost of data 
collection. 

In order to reduce the data collection cost, the lead 
time estimate can be computed in the same manner as for 
'Class I' items but for the class rather than individual 
medicines. This value can then be applied to each medicine 
in the group. 

To find the lead time estimate, a large number of 
latest orders should be selected and if the number of obser- 
vefions is sufficiently large, their distribution can be 
assumed normal on the basis of the contra! limit theorem. 
That value of lead time can be taken for all items which 
is 1 ange enough to include, say h-per cent, of the lead 
times of the items in 'Class-II' . This value can be readily 
estimated using normal distribution tables. The mean and 
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variance of 
on item 'i* 



a 



demend during leod time distribution for 
are given by the relationships. 


* 

Q. 




^m. 


' i 


■^i(h-percent) '* 


0 ^ 

if i(h-percent) 


(3.26) 


(3.27) 


ine order point can be found as before using Gamp and 
Mo'idel*s eztension of Tchobycheff * s inequality. As these 
items are important and their annual rupee usage is lower, 
we can opt for higher service levels. 

3.10 Control Procedure for Glass III Itemsi^ 


Those items have low overall! importance and hence the 
uss.ge of adaptive control limits does not seem to be economical 
as it involves considerable data collection effort. Therefore, 
fixed limit controls have been used for these medicines. Fo 


forocasting is needed on a period by period basis and th© 


lot sizes have been computed using the classical BOQ formula 
which states that for an item 'i*, 




(2 C^g) X 100 


^i 


%c 


(3.28) 
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wlier G , 


Q. 


°1S 

p. . 

1 

Pic 


= optimsl lot size, 

= azuiual demand in units, 

= set up cost per order 
= price per unit, and 

= carrying charges expressed as percentage of 


’Pi’ 


I'his lot siae is, of cours<.. , checked against the expiry time 
of the medicine. 

In calculating the reorder point, instead of concen- 
trating on accurate lead time estimates, it is more economical. 
to provide larger buffers. The buffer is provided for a 
period 'T-„' and is obtained by. adding this period to the 
*h-percent’ confidence limit lead time value as follows; 


R. 

10 


* 


0* 


’^i( la-percent) 




in 


is 


or 


R. 


10 


Q. 


m 


I- Pi(h-percent ) ^ ^is ^ 


( 3 . 29 ) 


where ‘T. ' is the fixed time factor whose value is assigned 

IS 

by the management. A greater value of ’T^^’ is used if the 
lead time for the medicine is long and variable and ii a 
stock-out of the medicine is particularly undesirable 
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froEi some other considerations. 

5 . 11 Summary 

In this chapter an inventory classification scheme 
based on fuzzy clustering concepts has been developed. 

Further a detailed an count of the inventory control proce- 
dures to be used for various classes of medicines has been 
presented. The proposed inventory classification and control 
procedures have been tested using a case study on medicine 
inventory i The fuzzy inventory classifica.tion scheme is 
compared with the classical ABC analysis. The data for the 
case study has been collected from tho Haaith Centre, Indian 
Institute of Technology, Eanpur. The results are presented 
in Chapter IV. 



CHiPIER IV 


A CASE STUDY s jRBSULTS, DISCUSSIONS AND THE SCOPE 

POR FURTHER WORE 

In order to validate the methodology developed during 
the course of this work, the medicine inventory control prob- 
lems of the Health Center, Indian Institute of Technology, 
Kanpur were studied. The discussion on hospital medicine 
inventory control problems presented in the previous 
cha.pters was also to an extent motivated from the observa- 
tions made there. 

The currently used medicine inventory control 
procedure at IIT Kanpur Health Center is the same as described 
in Section (3.6), The specific problems encounteiiRl while 
designing a medicine inventory control system have been 
discussed in Section (3.3). 

Data has been collected for 119 medicines which 
represent a cross-section of the medicine inventory at 
IIT Eajapun Health Center. These medicines have an annual 
usage of approximately Rs. 4,32,000. Even as sinning that 
the inventory set-up and ca,rrying costs add up to 10 per 
cent of the inventory investment, it is hoped that through 
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proper control procedures a saving of upto 25 percent, tliat 
is Rs. 10,000, can be realised annually on this small cross- 
section of the medicine inventoiy. The second important 
advantage will be the attainment of higher service levels 
which are of utmost important to service organisations like 
hospitals. 

In the following sections results based on this 
case study will be presented. First a traditional IBO 
analysis has been performed which is followed by FIGS (Fuzzy 
Inventory Classification Scheme) analysis. The results 
obtained from both the analyses have been compared. A 
computer program has been developed in FORTRAN IV to compute 
the parameters required for implementing the designed 
control system for various classes of medicines. The 
details regarding this program and other programs developed 
during this work are given in Appendix I. 

4.1 ABC Glassification 

The procedure to perform annual rupee usage based 

ABC classification, given in Section (5. 4»1) > 
computerised. The listing of the computer program is given 
in Appendix II. The results of this analysis are as follows: 

1. All the 119 medicines, taken toother, have an annual 

usage of Rs. 4,52,237/ -• 
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2. The first 20 percent of the medicines, i.e. 23 medi- 

cines belonging to the A-class have an annual rupee usage 
of Rs. 307,520 which is 71.15 percent of the total. 

ine middle 30 percent of the medicines, i.e., 

35 medicines belonging to the B-cl8.ss have an annual rupee 
usage 01 Rs, 1,03,735 which is 24 percent of the total, 

last 50 percent of the medicines, i, e. , 61 medicines 
belonging to the C-class have an annual rupee usage of 
Rs. 20,982 which is 4.85 percent of the total. The class to 
which a particular medicine belongs is given in Table 4.1. 

Figure 5 depicts the ABO curve based on this analysis, 

4 . 2 Fuzzy Inventory Glassification 

The same set of medicines has been subjected to 
FIGS analysis. A questionnaire given in Figure 6 has been 
used to collect the infoimation about the various feature 
vector components. The first step in FIGS is choosing the 
initial centroids. 

4.2.1 Initial Gentroids 

Various sets of initial centroids have been generated 
by using the Ball and Hall approach described in Section (3.5.4). 
Figure 7 shows the variation in the number of centroids 
obtained with respect to the Inter-centroidal distance. It 
indicates that as the required separation increases, the 
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Pigure 6: Questionnaire used for data collection. 


Questionnaire 

1. Name of the Item: 

2. Cost per unit: 

3. iinnual consumption: 

4. ilnnual rupee usage: 

5. Expiiy time in months: ' 

6. Maximum percent quantity 

discount available: 


7. Set up cost per order: 

8, Carrying cost /item /month: 


9. 


Item type: V. Important 
Life Saving 


j Important 

i 

. -t 



Ordinary 



10. Substitution-: Never 

Possible 


j Sometimes 

I 

.( 


/ilways 



11, Market 

Availability: Difficult f j 

* ! 


Sometimes 
not avai- 
lable 


Always 



12. Shortage 
history: 

13. Storing 
method: 

14. Storage 
space 


Frequently 

short 

Special 

car^ 

Occupies 
too much 



Sometimes | 
short 1 



Some 1 j 

t 

c are j | 

i 

1_ j 

Ave rage 

— <^1 


No stock f 
outs 1 


No care f 


Low 
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Fig. 6: Form for lead Time and Demand Data 


Collecti 

.on. 
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J> TIMS 



iemii:d history 

No. 

T1 

T2 

T3 


MONTH 

YEAR 

YEAR 

YEAR 
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JAN 




2 





FEB 




3 





MAR 
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, 


4 





APR 


.. 


5 





MAY 




6 



r II 11 II r. 


dUN 




7 





JUl 


.... 


8 





AUG 




9 





SBP 




10 





OCT 




11 





NOY 




12 





DBG 




13 





TOTAD 


1 : 

1. 

14 
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15 
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r* — 
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16 





17 


1 - .... . 
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18 


t. 

i 


19 
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L, 



20 


|- 
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! 



T1 = Time when shortage noticed to time when 
purchase requisition accepted. 


T2 = Time of acceptance of purchase requisition 
to time when order mailed. 

T3 = Time f rom mailing of orders* to items ready 
on shelf for distrihution. 

T * T1 + T2 + T3 
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number of centroids generated decreases. These initial centroids 
have been taken as input to the Puzzy ISODATA clustering 
algorithm. 

4*2.2 Partitioning bv Puzzy ISODATA 

Por a given set of ’c’ initial centroids, several 
fuzzy partitions have been generated by varying the values of 
the exponent ’m'. hqual weights have been assigned to al 1 
the feature vector components during this process. The 
experiments have been conducted for c = 2,3 and 4 and both 
the measures of partition quality, H and P , have been 
evaluated for each fuzzy c -partition. Table 4.2 lists values 
of Pj, , ^c partitions obtained for various 

values of 'c* and ’m', Figures (8), (9) and (10) depict the 
variation in P and H with respect to ^m’ for c= 2, 3 > and 4 
respectively. It can be seen from Table 4,2 that with 'c* 
fixed, — ^0 and P . — ^1 as 'm* — > 1 , This is expected since 

and the partitioniug is optimal for only when it 
lies in M , the hard partition space associated with X [22], 

Thus as the value of *m*’ increases the fuzziness of the 
^■*»partitions of X obtained with Puzzy ISODATA also increases, 

low scanning Table 4.2 for minimum which is a 
better measure of quality [22], we observe that the entropy 
is least for c = 3 for every value of 'm*. Prom this we 
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infer ■th.o.t ihe most plausible choice for primary substructure 
in X is 3 subclasses. As no theoretical result suggesting 
an optimal choice of *m’ has been developed to date [22], 
further discussion is based on, « = 3 and m = 1.25, which 
appear to be satisfactory for our purpose. 

The initial centroids for starting the Puzzy 
ISODATA process arei 

= ^0.75, 0.93, 1, 1, 1, 0.5, 0.5, l^ 

^2 = [ 0.12, 0.27, 0.45, 0.36, 0.41, 0.36, 0.44,0.52] 

v^ = {o.Ol, 0.75, 0, 0, 0, 0, 0, l] 

The final centroids which have been obtained at the end of 
the classification process are: 

= 1o.2842, 0.6196, 0.4768, 0.2433, 0.4619, 

0.3572, 0.5032, 0.8790] 

= I 0.0412, 0.1162, 0.6850, 0.6569, 0.6098, 

0.5595, 0.5997, 0.4379 ] 

V, = {^.0499, 0.1139, 0.0721, 0.0545, 0.0613, 

0.0589, 0.1449, 0^2405 } 

These centroids essentially represent the three subclass 
present in the data set. The first cen'jjroid has distinctly 
higher annual rupee usage and expiry t4^® components (i.e., 


ill 
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first two components of the feature vector) as compared to 
the other two centroids. Moreover, the values of its other 
components are- also substantial. This clearly indicates 
that the firat centroid represents the group of those 
medicines which have a higher annual rupee usage and shorter 
expiry duration as compared to the medicines belonging to 
the other two classes. The other components of this centroid 
also show its considerable ovorall importance. Therefore, 
we conclude that the medicines belonging to ’Glass-I’, 
represented by this centroid, are important from l consi- 
derations and should be controlled by the best procedure. 

The second centroid has low annual rupee usage and 
expiry time components but its other components advocate 
that substantial importance be assigned to the group it repre- 
sents. In fact, this centroid essentially represents the 
group of those medicines which have lower annual rupee usage 
and larger expiiy period but are otherwise quite important. 

A good number of injections belong to this group. However, 
there are some medicines in this group which have low 
expiry time. This is an example of misclassification. But 
such medicines can be detected as they have a lower me mb e3>- 
ship grade in this group as compared to the other medicines 
belonging to it. Two such examples are CAICIUH SAHBOZ 
INJEGTIOH and DEGABHOH liTJEGTIOH. Anyway as we are making a 



115 


check during the lot sizing ceaculations to see that an 
order is not placed for a period greater than the expiry- 
time of the medicine, the errors due to this misclassification 
will not he serious. As these medicines are important and 
ha-ve low annual rupee usage and high expiry time , we can 
opt for higher service levels for them without sacrificing 
much in terms of costs. 

The third centroid has distinctly low values of 
all the feature vector components and it oh-vLously represents 
the medicines which are unimportant from all considerations, 

A relatively loose centrol procedure can therefore be used 
for this group of medicines. 

We would like to emphasise here that any conversion 
of the fuzzy partitions to hard partitions should not be 
made at once by using some rule like the 'maximum membership 
grade rule'. Instead we should leave the solution fuzzy 
to identify the individuals in the data whose characteris- 
tics are shared by competing subclasses. When this has been 
thoughtfully considered, for all practical purposes the 
partitions may be converted to the hard c— partitions by 
using a suitable rule* This will greatly help in avoiding 
any misclassifications. 
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lie nav! oomrorted tile three fuzzy partitions In 
to hard partitions by using the aaxinui. membership grade 
rale and the resulting elassifioatlon has bec-n oompaKd with 
partitions obtained from the classical ABO analysis. 

4B0 Analysis v/.g iTTn.c; 


lysis. 


Table 4.2 shows the membership grades of the medicines 
in thethreo classes which have been identified by FIGS. It 
also shows the class ta which they belong to after the 
fuzzy partitions have been converted into haixJ partitions 
by applying the 'maximum membership grade rule' and the 
class to which they belong in a traditional 2 O-- 50-50 ASO 
classification.. Whereas in an ABC classification system the 
most important, important and unimportant items belong to 
the A, j 3 and 0 classes respectively, in FIG 3 the correspond- 
ing classes will be 'Glass I' , 'Class II' and 'Class IM* . 
Table 4.5 shows the percentage of medicines falling in each 
group of the two classifications. 

Table 4.5 « Comparison of ABC Analysis with FIGS 


Glassification Percentage of it ems in Class 
system , j 3 jjj or C 


FIGS 


52.77 


41.17 


'50 

26.06 
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The details regarding the changes in the membership of 
various items in the three classes are as follows j 

1. Out of the total of II9 medicines, no change has 
"takGn place in the classes of 61 medicines. 

2. No transfers have taken place from ’Glass 1’ to 
'Class II' and only two transfers have been recorded from 
'Glass A' to 'Class III', This implies that out of the 
total. 23 medicines belonging to 'Class A', 21 have retained 
their memberships in the most important class of FIGS. 

The two medicines transferred are ASGABIOl EMULSION and 
HSLOXON - 500 mg tablets. An examination of their feature 
vectors will reveal that they s.re really unimportant from 
considerations other than the annual rupee usage. The 
feature vectors are 

Ascabiol Emulsion = -[0.23, 0, 0, 0.5, 0, 0,5, 0, 1^ 

iledoxon-500 mg Tablet =[o.29, 0, 0, 0, 0, 0, 0, 1^ 

Hence these changes are quite justified and are according 
to the philosophy of selective inventory management. 

3. Sixteen medicines from 'Class B' have been promoted 
to 'Glass I*. These changes are also satisfactory since 
their feature vectors justify their claim for being 

controlled by better procedure. For example, Penzyme and 

Unienzyme Tablets which belongs to this class have a 
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feature vector, 

[.053, 1, 0, 0, .5, .5, .5, l} 

4. Pour medicines have been downgraded to ’Class III’ 
from Glass B’, Their overall low importance is evident 
from their featute vectors. The feature vectors are, 

LIV - 52 Tablet = {O.IO, 0, 0, 0, 0, 0, 0, I'J 

LIT - 52 Ibrops = { 0.09, 0, 0, 0, 0, 0, 0, 1^ 

Mejoral Tablet = [ 0.10, 0, 0 , 0, 0, 0 , 0, 1 } 

Sedonal Tablet = {0.05, 0, 0, 0,5, 0,5, 0, 0.5, o} 

5. Two medicines for which the feature vectors are given, 
below, have been promoted from ’Class C' to ’Class I’. 

Ajdexolin Liquid =[.018, 1, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5^ 
Glaxena Tablet = {,023,0.928, 0.5, 0.5,0,5,0.5,0.5,0.5| 

The feature vectors of these medicines also favour 
this change. 

6. Thirty four medicines have been promoted from ’Glass G’ 
to ’Class II’. A good number of them are injections which 
have low annual rupee usage but are important from other 
considerations. This is the peculiar chaise t eristic of 
'Class II ’ , and therefore, this change is also rational. 
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From the observe^tions above we conclude that the 
FICo has been yery successful in recognising the underlying 
importance groups in the medicine inventory and the changes 

brought about by .it in the tranditional mO classification 

* 

are justified. It has correctly classified the medicines 
on the basis of their overall importance in order to subject 
them to proper control procedures. 

4 • 4 Integrated Forecasting and Lot-Sizing 

A computer program, has been developed which forecasts 
the demands and computes the lot sizes and reorder levels 
for all the three classes of medicines as per the proce- 
dures described in Sections (3.8, 3.9, 3.10). It also plaints 
out rela.ted information on the investment required for 
each item and expected time for placing the next order. 

A sample output has been shown in Table 4.4. The program 
works on independent ordering policy for all medicines and 
whenever a replenishment order is to be placed, it needs 
to be run only for the required medicines. 

4.5 Computational Experience 

The greatest drawback of fuzzy clustering methods 
is a relatively large memo ly requirement. The core memory 
needed for FICS exceeded 16 £. However, with new and more 
powerful computer systems this . should not remain a handicap. 
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The ti'iio's teken by various progrsjns on IBM 7044/1401 system 
are as follows? 

(i) Claseicsl iiBC analysis — 32 seconds; 

( ii ) Initial centroids by Ball and Hall method - 46 seconds 
for all sets of centroids, 

(iii) j-'fizzy ISOBATii classification for a single value of 
’e’ and ’m' - on average 60 seconds, 

(iv) i'or lot size reorder point calculation - 30 seconds 
for 40 items. 

Though a, good amount of computer time will be spent 
in initially finding the number and types of classes present 
in inventory, this has to be done only when considerable, 
changes in da,ta profile have occurred. However, once the 
various subclasses have been established, the subsequent 
reclassification runs will not take much computer time. 

4*6 Scope for Further Work 

In the course of the development of the methodology 
presented iq this dissertation, it was felt that further 
research needs to be undertak.en in several directions. 

The potential areas of investigation and the possible 
extensions ares 
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1. The inventory control system is semi-computerisod 
in its present form. The records are still maintained 

in a ledger and manual updating and checking is required 
to ensure the timely replenishment , of medicines. As a 
further extension, the system can he fully computerised 
for automatic data processing. All records can he 
maintained on storage devices on computer and the transac- 
tions can he directly entered into the files for automatic 
updating. Whenever the stocks deplete helow the reorder 
point, exception reports may he printed out. 

2. It has been pointed out earlier in this thesis 
that the proposed inventory control system can he easily 
extended to suit the requirements of drug warehouses, 
medicine shops and retail stores. It will need only 
slight modification in the feature vectors. This type of 
research will help in further valida,tion of the proposed 
methodology. 

3. The applicability of more sophisticated techniques 
for forecasting the demands of a large number of items 
should he studied and if found suitable they can he 
incorporated in the proposed invoitorj^ control system. 

4. Joint ordering policies may he tried for medicines. 

Por better results, the dynamic lot sizing models, accounting 
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for varying rates of return and prices and price breaks 
can be incorporated into the proposed inventory control 
system, 

5. Extensile research needs to be xindertaken in the 
theory of fuzzy sots. As wo have pointed out ea,rlier in 
this chapter, so far no theoretical results are known to 
optimally choose initial centroids and the exponent *m' 
for fuzzy clustering. Further research regarding this will 
greatly reduce the computational burden and make the 
algorithms more effective. 

6, More attention should be paid to the selection of 
features because contrary to one's intuition, adding more 
features does not always lead to better performance of 
classification algorithms [52], ¥e have presented the 
results only for the case when all the features are assigned 
equal weights. Further research should be undertaken to 
study the weightage sensitivity of the algorithm. 

Lastly, the theory of fuzzy sets seems to have many 
potential applications in industrial engineering and attempts 
should therefore be made to use this approach in the areas 
where the situations are Intrinsically fuzzy. 



125 


4.7 Summary 

¥e havt. seen that the Fuzzy Inventory Glassification 
Scheme is more general and compatible with the concept of 
Selective Inventory Management. In this sense, it is 
much more superior than the classical ABC analysis. Further, 
with the introduction of scientific inventory management in 
hospitals, significant improvements can be expected in 
the present cost structure and service level. The theory 
of fuzzy sots itself shall prove to be a powerful tool 
in analysing the complex systems and it undoubtedly has 
great potential for practical applications in industrial 
engineering. 
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iiPPBNDIX I 


USER'S MANUia, 

In "this appendix the details of the various input 
and output parameters for the computer programs developed 
during the course of the work presented in this disserta- 
tion are given, Pour programs have heen developed in 
PORIiltiN IV and run on the IBM 7(^^/1401 Computer System at 
III Kanpur. The first program deals with classical ABC 
analysis and the remaining three programs correspond to the 
inventoiy classification, and cpntrol system designed in 
this work. Though the latter throe programs can be easily 
combined into a single program* <iue to available memory 
restrictions they have been retained as separate programs. 
The description of various FORTRAN variable names used in 
each program is gi^en at the beginning of the program; 
listing. Comment cards have been used liberally throughout 
the program to make it more readable. 

1. Classical ABC ilna lYSis Program 

FUNCTION? This program performs the classical ^IBG analysis, 
(i) The percentage of items to be put in each 
class has to be specified in a BATA statement in 
the beginning of the program. Olie percentage of 


INPUTS: 
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OUTPUTS 


items in class B aid C is denoted by /iPER, 

BPER and CpER respectively, 

(ii) The data cards are designed in such a way that 
the same cards can be used for the first three 
programs. The details of the data cards are given 
in the next section. Prom each data card the name 
of the medicine and its annual rupee usage are read 
as follows; 

Variable lames ; RilME (I, J), USAGE (I). 

Pormati 40A1, P5.0 

(iii) In its present form the program can handle 150 
items, but this number can be increased by simply 
changing the DIMBUSIOU statement. The program 
automatically computes the number of items for 
which the data has been fed. The end of the data 
file is indicated by the last card which has ’ * * 
punched in its 1st and 40th columns. 

The following information is obtained as output; 

(i) The input medicine inventory data. 

(ii) Results of the ABO analysis in a concise form. 

(iii) Name and annual rupee usage of items belonging 
to each class. 

/ ;r-r . -- 1.’; "• : ■ v'^ . 
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FUNCTION : It locates initial centroids ly Ball and Hall’s 
method for given values of the weighing exponent 
’EMM’ and weightages '¥( J)’ assigned to feature 
vector components » for several specified values 
of the inter centroidal distance ' DIST’ . 

INPUTS : (i) The value of the exponent ’EMM’ is given in 

the beginning of the program in a DATA statement. 
Its value can be changed by simply modifying 
this card. 

(ii) The first data card contains the total 
number of items, N, and the number of components 
in each feature vector, M. 

Variable Names: N, M 

Eoimat : 213 

(iii) The second data card contains the values 
of weights to be assigned to each component of 
the feature vector. 

Variable Name J W( J) 

Format s 8F5.2 

(iv) The next set of <3 at a cards contains the 
names and feature vectors for various medicines . 
G or^sponding to each medicine . there is one 
data card which is punched in the following 

format. 
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OUTPUT J 


Variable Fame: X(I, J) 

Format ; 40AI, F5.0, F3.0, 6F3.1 

The first 40 columns af the data card contain 
the name of the medicine, the next 5 columns 
it’s actual amual rupee . usage and the 
subsequent 3 columns the actual expiry time 
of the medicine in months. The rest of the 


columns have the values of the remaining six 


feature vector components punched as either 
0, 0.5 or IjD. However, this program does not 
read the. names of the medicines. These data 
cards are common for the first three programs 


and in each program only the relevant information 
is read. 

(v) iifter the data cards for the medicines, each 
of the subsequent 20 data cards contains one 
value of the intercentroidal distance ’BIST’, 


Variable Name: BIST 
Poimat F5,2’ 

file obtained from the 

computer output: 

bk 

distances and the corres- 
for them. 
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(vi) The character ‘ * » is punched on the 1st and 40th 
columns of the last datacard. This acts as a 
signal to the computer that the data file has been 
completely read. 


OUTPUT : The following information is printed out: 


(i) . Input data. 

(ii) Total number of medicines considered for analysis, 

(iii) lformalised feature vectors of all the medicines, 

(iv) The total number of iterations required for 
convergence, final centroids, names of the items 
and their membership grades in each fuzzy cluster 
and the values of partition quality measures, 

1-Pq and Pj,. 

4. Program for Combined Forecasting. Lot Sizing and Reorder 
1*0 in t Set t ing — 


FUNCTION jThis is a common program for calculating lot sizes 
and reorder points for the medicines belonging to 
the three classes. There is no restriction on the 
number of medicines in input data as the program 
works on independent ordering strategies, 

INPUT: (i) The values of K, N, ALPHA, ETRATE and SETUP are 

specified in the first DATA statement in the beginning 

of the program. 
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(ii) The values of Tl, T2, YARTl and 7ART2 are 
specified in another DATA statement. These values 
will be used for'Olass I' and 'Glass II' medicines. 

The following data cards are required for each 
medic inesj 

(iii) The first data card contains the following 
information. 

Variable Name: NAME, PRICE, EITIMB, ICLASS, FCURNT 
Format: 40 ;a, E4.2, E3.0, 213- 

(iv) For the medicines belonging to'Class I' and 
'Class II,' the rest of the Data cards are as follows} 

(a) Second data card: Variable Name: T3MBAN, VART3 

Format: 2F6.0 

(b) Next three data cards contain monthly demand 
data of last three years. Bach data card has the 
information for one year. 

Variable Name: D(E1, J) 

Format : 12 P6.0 

(c) The last data card contains current year's 
demand data for 'NCURNT' months in the same format 
as specified in (b). 

( v) The service levels for the'Class— I' and 'Class 127 
medicines are specified by two’ DOGIC AL IF’ statements. 

(vi) For the medicines belonging to'Class III J only 
a single data card, after the first card specified 
in instruction (iii) is required: 



ouTPag? ; 


“Variable lame: AMUKL, TTOTilL, SUPPER 
Pormat : 3P6.0 

(vii) The data cards for all the medicines are 
prepared as above, Afe a time data for a single 
medicine is read and processed by the program, 
(viii) The end of the data file is indicated by 
a card whose 1st and 40th columns have the 
character ' * * punched on them. 

The output contains the follov/ing information; 

(i) An exception report, if the current year has 
ended, asking the user to update the data file,. 

(ii) A report is printed for every item containing 

its name, optimal order quantity^, reorder level, 
unit price, cost of the lot in Rs,, number of 

month in which the next order is likely to be 
placed. The month for the next order is indicated 
by a number between 1 and 12. 

The FORTRAN listings of these programs are given 


in Appendix II, 
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APPENDIX II 


PRO GRAM IISTINGS 




IBFTC KMN 


CLASSICAL ANNUAL RUPEE USAGE A-B~G ANALYSIS PROGRAM 


NAME(I,J) • name of THE ITEM III 

USAGE! I } « ANNUAL RUPEE USAGE OF THE ITEM (I) 

N * NUMBER OF TYPES OF ITEMS IN INVENTORY 

NA » NUMBER OF ITEMS IN A-CLASS 

NB « NUMBER OF ITEMS IN B-CLASS 

NC - NUMBER OF ITEMS IN C-CLASS 

APER « percentage OF ITEMS BELONGINGING TO A- i 

BPER » PERCENTAGE OF ITEMS BELONGING TO B-CLASS 

CPFR * PERCENTAGE OF ITEMS BELONGING TO C~CLASS 


DIMENSION NAME tlSOf 40), USAGE! 1501 


DATA APER,BPER,CPeR/20.,30,,50«/ 
FORMAT!/, 5X,l2l!lH-)) 


DATA ISTAR/IH*/ 

READ 1004, (NAME!I,J),>1*40),USAGE!I> 

IF*^n^NAME(iaUEQ«lSTARJ .AND. (NAME! 1, 40y.EQ,ISTAR) > GO TO 2 

I » I + I ' ' ' ' ' ' ' 

GO TO I 


PRINT 1000 

FORMAT!/, SXfUKlH*)) 

PRINT 1006 

FORMAT t7t42X,*THE INPUT 
PRINT 1000 
PRINT loot 

F0RMAT!/#25X,*ITEH NO.*, 

I SAGE*) 

PRINT lOOO 
DO 3 ' I 

PRINT 1008 , 1 # !NAME!I, 4), 


MEDICINE INVENTORY DATA IS AS FOLLOWS*) 


♦NAME 3F THE ITEM*, 30X, ♦ANNUAL RUPEE U 





o o o o 



CONTI NUE 

PRINT 1000 


arranging the items in descending order of usage 


NN * N - I 
ISW •* 0 

DO 130 I * 1,NN 
IF CUSAGECn.GE.USAGEn+ll I 
ISW « ISW ♦ 1 
TFMP « USAGE! I) 

USAGE (U « USAGE! I+l) 

USAGE !I+IJ « TEMP 
DO 120 J » 1,40 
TEMP » NAMEilfJ) 

NAME! I, J) * NAME(I+1,JI 
NAME! 1+1,4) » TEMP 
CONTINUE 

CONTINUE v; : 

IF !ISW*Ne.O) GO TO 105 


DIVIDING THE ITEMS IN TO A,B, AND C CLASSES 


NA ■ IFIX! (FtOATCN)*APER)/lOO. ) 
NB « IFiXf !FLQATtN)TBPER)/lOO, ) 
NC » N - ( NA + NB ) 

ASUM • 0.0 

DO 201 I * IfNA 

ASUM * ASUM + USAGE! I) 

CONTINUE 

BSUM .0.0 '■ ,,, 

DO 202 I « 1#NB 

INA * I + NA' :Si,- ■ 

BSUM * BSUM + USAGE! INA) 

CONTINUE # '4 ■ 

CSUM . 0.0 ; '' ''§ m 

DO 203 I • IfNC > ' ‘.'.0 
lAB . I + NA ♦ NB 
CSUM . CSUM USAGETIAB) 

CONTINUE " ' 

TSUM . ASUM > BSUM + CSUM : 

PH inY 1010 ' . ‘ ' 

FORMAT !/+40X»*THE^ RESULTS OF;-^^^BH 

AS PER i”? A SUH AIOO. ) /TSOH-Wi; • 

- I RSUM^lOO* ) /TSUM ; 

CSPER - 100.00 - (ASPER ♦ 6SPER ) 

PRINT 1011, APER, ASPER 


1011 

1012 

1013 

1014 


301 

1017 

1010 


PRIMT 1012tBPE«tBSPER . ■ ■ 

PRTfjT 1013,CPER*CSPER 

FORMAT I /,25KfF6*2t2Xt*PERCENT ITEMS BELONGING TO A-CLASS HAVE ANN 
lUAl RUPEE ISAGE OF*f2XtF6.2t2X, ♦PERCENT*) 

FnRPATC/,25X,F6*2f2Xf*PERCENT ITEMS BELONGING TO 8-CLASS HAVE ANN 
lUAL RUPEE USAGE OF*t2X,F6*2,2Xf *PERC£NT*) 

FORPAT(/f25X,F6.2,2X,*PERCENT ITEMS BELONGING TO C-CLASS HAVE ANN 
lUAL RUPEE USAGE OF*t2X»F6*2,2X»*PERCENT*) 

PRINT 1000 
PRINT 1014 

FORMAT (/»40X»*THE FOLLOWING ITEMS BELONG TO A-CLASS*) 

PRINT 1000 
PRINT 1007 
PRINT 1000 
DO 301 I*lfNA 

PRINT lOOSflf lNAMElI*J)tJ*lt40l,USAGE(I) 

CONTINUE 
PRINT 1001 

PRINT lOlTfASUR 

format t/»40X**TOTAL ANNUAL RUPEE USAGE OF THIS CLASS OF ITEMS **» 
I lXt*RS**tF10.2) 

PRINT 1000 

PRIHT 1018 

FORMAT C/f40X»*THE FOLLOWING ITEMS BELONG TO 8-CLASS*) 

PRINT 1000 
PRINT lOOT 
PRINT 1000 
NAl • NA ♦ 1 
NANS • NA ♦ NB 
00 302 1 • NAltNANB 

PRINT* lOOBfIBf (NAMEdf J)»J*1»A0)*JSAGE(I) 

302 CONTINUE 
PRINT 1001 
PRINT 1017, BSUF 
PRINT 1000 

PRINT 1019 .nu»Mc TTFM<; RFL0N6 TO C-CLASS*) 

FORMAT!/, 40X, ♦THE FOLLOWING ITEMS BELUNb 

PRINT 1000 
PRINT 1007 
PRINT 1<J00,:" ' 

NANBl * NA N8.*' 1' , 

OO 303 I * NANBl, N 

IC ■ I * n. i»1.40)tUSAGE{I) 

PRINT 10G8,IC,lNAMEtI, •))*'» 1» 

303 CONTINUE 
PRINT 

PRINT lOi r,^ 
print IIJOO 

,PRINT'MEZQ*tSUM 


1019 


1020 


I annual rupee usage of all ITEMS **,1X,*RS.*, 

PRINT 1000 

STOP ^ , , 

END 



oooo h- f-* ^£) I-* ooooooooooooooooooonoooooon 


IBFTC ^'AIN 

THIS PROGRAM LOCATES INITIAL CENTRIOOS BY BALL AND HALL METHOD 


NOTATIONS USED 


N = NUMBER OF TYPES OF ITEMS IN INVENTORY 
M = NUMBER OF COMPONENTS IN ANY FEATURE VECTOR 
X(I,J) * COMPONENT (J) OF FEATURE VECTOR OF ITEM 11) 

X(I,1) = ANNUAL RUPEE USAGE COMPONENT 
Xflf2) = EXPIRY TIME COMPONENT 
XCI»3) = LIFE SAVING IMPORTANCE COMPONENT 
X{I,4) » SUBSTITUTABILITY COMPONENT 
X{I,5) = MARKET AVAILABILITY COMPONENT 
X(I,6) = STOCKOUT HISTORY COMPONENT 
X(I,7) = STORAGE SPACE COMPONENT ’ 

X(I,8) » STORING METHOD COMPONENT 

V(K,J) = COMPONENT (J) OF CENTROID VECTOR OF CLASS CK) 

W(J) = WEIGHT AGE ASSIGNED TO COMPONENT (J) OF F,V. OF ITEM 
EMM = WEIGHTING EXPONENT 

DIST = SPECIFIED I NTER-CENTROI DAL DISTANCE 


DIMENSION XI150,12)*V(60,12)»D<30) *Wf 12) ,VAR ( 12) 
DATA EMM/1.25/ 

READC5,1001) N,M 
001 F0RMAT(2I3) 

WRITE(6,1001) N,M 
READiBiDDR) (W(J),J»lfM) 

99 FORMAT! 8F5. 2) 

WRITE(6t999) (W(J)tJ=lfM) 

READ! 5f 1009) ( !X<I,J) t J*lf M)f I=1,N ) 

009 F0RMATl40X,F5.0tF3.0t6F3.1) 

WRITE(6,1109) !!X(I,J)f J=ltM), I = lf N) 

109 F0RMAT(40XtF8.0t7F6. 1) 

EPP * 1.0/(EMM-1.0) 

ADJUSTMENT OF EXPIRY TIME IN FEATURE VECTORS 


GREAT * 0.0 
DO 4 I = IfN 

IF (GREAT. LT.X(If2))GREAT * X(It2) 

CONTINUE 

DO 5 I = IfN 


4 



X(I,2) = ( GREAT - X<I,2) ) / GREAT 
5 CONTINUE 

C 

C normalization of feature VECTOR COMPONENTS 

C ... 

c 

DO 35 J*1,M 

IF ( J ,EQ. 2 ) GO TO 35 
GREAT = 0,0 
DO 25 1=1, N 

IF (X{I,J) .GT. GREAT ) GREAT = XtI,J} 

25 CONTINUE 

DO 30 1=1, N 

X(I,J) = X(I,J) / GREAT 
30 CONTINUE 

35 CONTINUE 

C 

C THIS PORTION LOCATES INITIAL CENTROIDS 

C IT USES BALL AND HALL METHOD 

C — — 

C 

CALL VARIAN(X,N,M, VAR) 

DO 300 III = 1,20 
REA0(5,1021) DIST 
1021 F0RMAT(F5.2) 

WRITE(6,1150) DIST 

1150 FORMAT!//, 5X,*F0LL0WING ARE THE CENTROIDS FOR DIST =*, F8.3) 

DO 210 J=1,M 
SUM = 0.0 
DO 200 1=1, N 
SUM = SUM + Xtl, J) 

200 CONTINUE 

V(l,J) = SUM/FLOAT(N) 

210 CONTINUE 

K » 1 

DO 265 1=1, N 
DO 240 JK=1,K 
SUM = 0.0 
DO 230 J=1,M 

SUM = SUM + W(J) * ! ( IXCl, J)-V( JK, J) )**2) / VARCJ) ) 

230 CONTINUE 

IF ( SUM .LE. l.OE-06 ) GO TO 265 
DtJK) = (SUM) ♦♦ EPP 
240 CONTINUE 

DO 250 JK=1,K 

IF ( DtJK) .LT. DIST ) GO TO 265 
250 CONTINUE 

K = K + 1 
DO 260 J = 1,M 
V{K,J) = X{I,J) 



260 CONTINUE 

265 CONTINUE 

WRITE (6,1170) ( lV(KJ,jy,J=l,M),KJ=l,K) 
1170 FORMAT!/, 20X,8F8«5) 

300 CONTINUE 

STOP 
END 

IBFTC VARIAN 
C 

SUBROUTINE VARIANIX, N,H,VAR) 

C 

DIMENSION X(150,12),VAR(12) 

C 

DO 320 J = 1,M 

SUMl = 0.0 

SUM2 = 0.0 

DO 310 I = 1,N 

SUMl * SUMl + X{I,J) 2 

SUM2 = SUM2 + XM, J) 

310 CONTINUE 

SUMl = SUM1/FLCAT(N) 

SUM2 = (SUM2/FL0AT(N) ) **2 
VAR(J) » SUMl - SUM2 
320 CONTINUE 

C NORMALIZATION OF VAR(J) 

GREAT = 0.0 
DO 330 J » 1,M 

IF ( GREAT .LT. VARlJJ ) GREAT = VAR(J) 
330 CONTINUE 

DO 340 J * l#M 
VAR(J) = VAR(J) / GREAT 
340 CONTINUE 

RETURN 
END 
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IKVHNTCRY 


CLASSIFICATION SCHFME 


I'rcTr,^; j-p INVENTORY CCNTRQL SYSTEM 
SYSTE!^ IS eased ON FUZZY SETS THEORY 
•'UZZY zsCDATA IS HPJNG USED FOR CLUSTERING 

-fptc prcgpAw classifies the items for effective ccntrol 


K r T A T I C ^ S USED 


K = 

M = 

K = 

xn? 

X( I, 

X( I f 

X{ I, 

X(T, 
X{T , 
X( I, 
X{ I, 
X{ I, 
V( K, 
V\C J) 
Ud, 
VI {J 
NAME 
FPSI 

FMM 
FC = 
1 1 « •“ 
HC = 
VAR{ 


MUMPER OF TYPES OF ITEMS IN INVENTORY 
NUMBER OF CGMPCNENTS !N ANY FEATURE VECTOR 
^^UMP.ER OF CLASSES REQUIRED TO BE FORMED 


( I ) 


J) = CCM'^OKENT (J) OF FEATURE VECTOR OF ITEM 

I) = annual rupee usage COMPONENT 

?) = EXPIRY TIME COMPONENT 

3) = LIFE SAVING IMPORTANCE COMPONENT 

4) = SIBSTITUTABILTTY CCMPONENT 

5) = MARKET AVAILABILITY COMPONENT 

6) = STOCKCUT HISTORY COMPONENT 

7) = STORAGE SPACE COMPONENT 

8) = STORING METHOD COMPONENT 

J) = COMPONENT (J) OF CENTROID VECTOR OF CLASS 
= We .TGHTAGE ASSIGNED TO COMPONENT ( J ) OF F.V 

K) = MEMBERSHIP GRADE OF ITEM{ 1) IN CLASS (K) 

,J). = THE SET OF NEW CENTROIDS REPRESENTING CLUSTERS 
< I , J) = NAME CF THE ITEM I 1 1 

L = CCNVERG*=NCe CRITERION FOR PROPER CLUSTERING 
= WEIGHTING EXPONENT 
PARTITION COEFFICIENT 
FC) = LOWER eCUNC OF ( FC ) 

NORMALIZED ENTROPY 

4) = VARIANCE OF THE (J)TH, F.V. COMPONENTS 


(K) 

OF 


ITEM 


DT 'MENS! ON X { 150, 1 2 ) » V ( 3, 1 2 ) , W ( 12 ), U( 1 50, 3 ) , PT( 12) ,GRDt 
CTM^MSTCN Vl{3,!2),SUNUR(150,12),SUPER(12),TUt3) 
DTMFNSTCN NAME f 150,40) , VAR ( 12) 


3) 


1 ^ 


PUT DATA 


ICOO 



n 'D n n 


1 r r 

1 fpp 

] r 0*’ 

r 

1 

IC'^A 

2 

ICO*^ 

10^6 

1G07 

1C08 

3 

C 

A 


r'T^ ^/n>/l,25/ 

IH*) ) 

inri » v,Kt EPSR 
FrF’'"\Tf 2TA fFR, 5) 

RO?.? {(VCKUJ), J=UP),K1=1,K) 

FCPf'Mf /.FF.f, ) 

P-EAr 1^03, {W{ J), 

ppp AT ( n F5 ,2 ) 

RFUT OF FEATURE VECTORS 

T = ! 

CATA TSX/^r/IH*/ 

PFAP lOCAj {AANP(I,J),J=l,AO),{X(ir,Kl),Kl=l,M> 
FCR‘^AT{40AljF5,OtF3.0,6F3, 1) 

TF { riAMFCT jl) ,EG« ISTAR) .AND. (NAME{I, AO) .EO. TSTAR))GO TO 2 
! = ! + 1 
GC TO 1 
N = T - 1 

PRTNT 1005 EPSTL 

FnRfMT{ /,7X,*NC. CF TYPES OF ITEMS IN INVENTORY =« , I A, 7X, =«‘NC. OF 
FEATURE VECTOR CCMPONENTS ==^, 1 4, // ,7X, *N0. OF GROUPS FORMED =* , 
I4,?5X,*CCNVFRGENCE CRITERION EPSIL ==«',FB.5> 

PR TNT lOOG 
PRINT lOOfc 

FCPMAT(/,42X,’!'THE INPUT MEDICINE INVENTORY DATA IS AS FOLLOWS*) 
PRINT inoo 
PRINT 1007 

FORWATI/f 5X,*ITEM MO. * 1 7X , *NAME OF THE I T EM*, 5 OX , *FEATLR E VECTOR 
*) 

pRT^r 1000 
CO 3 1=1, N 

PRINT 1008,1, (NAME(I,J) ,0=1,40), {X{I,J),J = 1,M) 

FORMAT! €X,T4,10X,40A1,5X,F8.0,7F6. 1) 

C0N""INUE : 

PRINT ICOC 

ADJUSTMENT CF EXPIRY TIME IN THE FEATURE VECTORS 

great = 0,0 

DO 4 1=1, N 

!F (GREAT .LT. X{I,2) ) GREAT = X{1,2) 

CONTINUE 
DO 5 1=1, N 

X(I,2) = (GREAT - X(I,?) ) / GREAT 
CONTINUE ' _ 


NOR'^ALJ ZATTCN CF THE WEIGHTAGE VECTOR 


GRCftT - fi,G 

DO 10 J =1,M 

f W{J) .GT. GREAT ) GREAT = WIJ) 



1C 

cr 1 ■ j=i , 
v ? J) W ? J ) / GRE/^T 
,r CnNTT'JlJE 

PilTM 101 1 

icn 55Xs^THE WEIGHTIiGES ARE AS FOLLOWS*) 

lC12jfW( J),J=1,M) 

1C12 FrRVAT{30X,«FR.'^) 

Pn’M ICGC 

c 

C PR’-RTTOG the tnitIAL centroids 

PR -NT 1013 

1C13 FCP^'AT{/,E5X,*THF INITIAL CENTROIDS ARE AS FOLLOWS*) 

PR^Nf^ 1000 

1014, f {V(K1, J), J=:1,M),K1 = 1,K ) 

1C14 FrRwAT(20X,8F12.8) , 

PRINT KjCC 
C 

r 

C NPRHALT ZATTON CF THE FEATURE VECTOR COMPONENTS 

C — ' — 

c 

CO O'" J=1,M 

TF (J .EC. 2) GO TO 35 
GREAT = 0, C 

no ?5 T=1,N 

IF {X(T,J) .GT. GREAT ) GREAT = X(I,J) 

25 CONTINUE 

CC 30 1 = 1, N 

X{T,J) = XfIjJ) / GREAT 
3C CONTINUE 

35 CONTINUE 

C 
C 

PRINT 102C 

lC?n FCR^'AT{/,45X,*THE FEATURE VECTORS ARE AS FOLLOWS* ) 

PRINT 1000 
PRINT 1021 

1C3I FORMAT! /,5X,*ITFM NUMBER*, AOX,*FEATURE VECTOR COMPONENTS*) 
PRINT ICOC 
CO 50 1 = 1, N 
CO AO .J=1,M 
PT{ J) = X{ I, J ) 

AC CONTT mjE 

PP’KT 1022,1, (PT(JJ),JJ=l,M) 

1C22 FORMAT! ’’ X , TP , 7X, 6F 1 2 . 5 ) 

" r CO KTI MJE 

C 



c 

c 


c 


62 


ft A 

68 


6 *7 


7 0 

71 

72 


73 


SC 

C 

c 

c 

c 


CCf^PUTTNG THG ^^E^BERSHIP GRADES 


C&LL VARI A^{X^^,MfVAR) 

FPP - 1,0/ (EWf^-l.O) 
i^CtNT = 0 
Cn\■^T^!UE 

¥Clf^'T = KCUNT + 1 
rC '’0 ^ = 1 ,N 
DC ^'2 KK=1 tK 
JGT = kK 
T=ST = C.C 
CC J = l,^' 

-TK-r-r ^ A^S (X{ I ,J)- V( KK, J) ) 4 test 

■’Fd^ST.LE. 1. 00-15 ) GO TC ’72 

fOMTMJE 

rn 60 KK=1,K 

GRffKK) = C.O 

nr (A 

GRCfKK) = GRDTKK) + W C J )=«'{{ C X (I , J ) -V { KK, J )) ) /VAR { J ) ) 

OrK-TKitF 

Gr>n(KK) = «GRD(KK)) EPP 

OF ( GRC(KK) ,IE. l.OE-20 ) GO TO 71 

CCMTNUE 

SUN = 0,0 

CO iFO kk = i,k 

SUN = SUN 4 1,/GRC(KK) 

continue 
00 ”0 kk=i,k 

Ufl,KK) = 1.0 / (GRD(KK)*SUM) 

CONTINUE 
GC TO 80 
JOT = KK 
CO 72 KK=1,K 
L ( I,KK) = 0,0 
CON’^TNUS 
iJ(!,jOT) = 1»0 
CCKTINUE 

TH?S PORTION CALCULATES NEW CENTROIDS 


CC 70 KK*1 ,K 
COKP =0,0 
Cn 0-4 T = 1,K 

CDAR = CCNR 4 {IJ ( I ,KK ) )=5^*ENM 

ccm:^ Nuc 

CO 86 ■' = l,N 

cn j=i , N 


3 A 





SU\URf^7J) = 

8 5 CCMTN'UE 

86 CON”' ^illF;■ 

HQ ■ J=l, K 
SUP -■>, t J) = 0.0 
CC 87 T=1 jK 

'•nP'EPfJ) = SUNljRfJ,J5 + SUPER(J) 
qT C'^NTIUUE 

Cn\TINIJE 
CC 80 J = 1,1^ 

VI ?KK j, J) = SUPER( J)/CDHR 

Q q CC6 TT ^^l IE 

9C crc’^'^ \UE 

c 

C TES'^ FOR CC''JVERGENCE 

c — 

c 

CO 100 KK=1,K 
CON! = O.C 

no 91 j=i,N 

CGM = (VlfKK,J) - V(KKfJ))^«2 + CCNl 
01 CCK’^TN'SE 

CON = SORT (CONI) 

TOfCON .GE. 5PSTL) GO TO 110 
ICO CCNTTMJE 

GO ■'0 120 

110 CO 115 KK=I,K 

CO 114 J=1,M 
VlKKjJ) = VlCKKjJ) 

114 coN'^T'^ne 

l 15 CONTINUE 

GO TO *"5 ' 

c 

C “'O^TRUCTICNS for PRINTING 

C — 

c 

120 print 1 100 

1100 format? //, sox, ’^ITERATIONS CONVERGED - ANALYSIS SUCCESSFUL*) 
PRINT 1000 
PRINT 1120 , KCUNT 

1120 FORMAT!//, 2SX,*T0TAL NUMBFR OF ITERATIONS REQU IRED =*, 1 3 ) 

PRINT 1000 

FR NT 1 1 2 1 

1121 FCRM/^t? /,*5CX,*THE FINAL CENTROIDS ARE AS FOLLOWS*) 

FKTN'T ICOC 

PR’^NT U25,( (V1?KK,JJ), JJ=1,M),KK=1,K) 

11?-; FCRVAT(/,8F12.8) 

PR’'NT ICOC . 

rjD T K'r- I ] ■^C ■ ' 



1130 FCR^^AT?//^^OX,>^M^^^egRSHIP GRADES ARE AS FOLLOWS*) 

ICOC 
PrrL,T 1 ) 32 

113? FORMAT? //,5Xt*ITEF NUr^BFR*,l£X,*ITEM NAP^E*, 30K,*GR0UP T*, 

1 1:5X,*GQCUF T T * , 1 S X , *GROUP III*) 

PHT^!T icCG 
CO 130 T = 1,N 

CC 123 KK= 1 ,K : 

-i.K'-'K) r. IH!,KK) : 

12 = 

114C ,I , (N6»^E{I, J) , J= 1,40 ), (TUCKJ ),KJ = ltK ) 

1141 FCRVAT»/,7X,T5,PX,40A1,3F21,8) 

I 3n rnN'^'^NME 

PR'^''!T icnc ^ 

cm C:U^LTMU,N,K,FC,HC) 

FCL = 1«'^ - FC 

ppff^;- 11SC,EMM,K, FC, FCLjFC 

II = 0 FDRWMf //, lOX, ^WEIGHTING EXPON E^^Tl EMM ) =*, FB.4, 10X,*M0. OF CLUST 

IRS(K) =*,T3 ,10X ,*P^RTITICN COEFFICI ENT{ FC ) =* , F 8, = , / /, 1 CX, *LOWF R ' 
70UKC CN (HC) = (l.~FC) =*,F8.5 , lOX , *NORM AL IZ ED ENTROPY f HC ) =*, 

3 F3«C5 I 

PR” NT ICQC 

S T C P : 

ENC ^ 

C 
C 

IPFTC VART /\N 1 

C 

SUBRnU'””NE VARlANi{X,N,M,VAR) ; 

CTMENSTCN X(150,12),VARa2) 

CC 320 J = 1,M 

SU^l = 0,0 

SIJVP = O.G 

CC 310 ! = 1,N 

SUMl = SUM + XCIfJ) ** 2 

SUM2 = SUV? + X{ I, J) 

310 CONTINUE 

SUVl = SUVl/FLCAT (N) 

SUM? = <SUM2/FLCAT{N) ) **2 
VARf J) = SUMl - SU«2 
320 CONTINUE ■ 

C NCR^ALIZATICN CF VAR{J) 

GREAT = 0,0 
CO 330 J = 1, M 

TF ( GREAT ,LT. VAR{J) ) GREAT = VAR(J) 

330 CCN”T MJE 

CC 340 J = 1,M 
VAP'J) = VARIJ) / GREAT 
0 CONTI NUE 



TMP 


c : 

IRF"r. CU.AI.T! 

C 

c 

^IJPUOUTTNF QUaiTnUfKtK, FC,HC) 
PT wir.,!C7r.N UC130,3) 

FC = n.o 

rr 

T ■ = n , r 

rn ■'o KK = I,K 
‘ Tc'^'P = l; ( T , KK ) ^*2 + TEMP 

1C CC^■^■r^||E 

FC FC + TEMP 
F COM'^TNUE 

FC ^ FC/FLCAT{K) 

HC = 0,0 

nC 20 KK = 1,K 

CO 20 ^ = IjN 

?e HC = HC + U{T ,KK)«flLCGnJ{I»KK) ) 

FC = -fFC/FLCATfN) ) 

pcTijriv 

CNC 



MAIN 


PROGRAM FOR FORECASTING - DYNAMIC LOTSIZING - REORDER POINT CALCULftTIONS 


THIS PROGRAM CALCULATES VARIOUS PARAMETERS REQUIRED FOR 
MEDICINE INVENTORY CONTROL FOR ALL CLASSES OF ITEMS 


K = NUMBER OF YEARS FOR WHICH PAST DEMAND DATA IS AVAILABLE 
N = NUMBER OF PERIODS IN A YEAR 
ALPHA = SMOOTHING CONSTANT 

AVERGEtJ) = MONTHLY AVERAGE OF (K) YEARS FOR THE PERIOD CJ) 
YEARAVd) = AVERAGE MONTHLY DEMAND IN YEAR Cl) 

RATIO(J) = DEMAND IN LAST MOST YEAR IN PERIOD! J) /BASE VALUE OFCJ) 
RBAR = SMOOTHED VALUE OF THE RATIO CJ) SERIES 
TREND = TREND PER MONTH 

BASE SERIES CONTAINS THE SEASONALITY INFORMATION 
RATIO CJ) SERIES INFORMS ABOUT SHIFT FROM BASE CJ) 

CRBAR) REPRESENTS THE SMOOTHED RATIO VALUES 
TREND GIVES THE TREND-CORRECTIONS FOR COMING YEAR 
DCI,J) » DEMAND FOR PERIOD CJ) OF YEAR Cl) 

FRCST(J) = FORECAST FOR MONTH CJ) 

NCURNT = NUMBER OF PERIODS OF CURRENT YEAR FOR WHICH DATA AVAILABLE 

PRICE = PURCHASE PRICE OF ITEM IN RS./UNIT 

SETUP = SET UP COST IN RS. /ORDER 

CARCST » CARRYING COST IN RS. /ITEM/PERIOD 

RERATE = PERCENT RATE OF RETURN ON CAPITAL FOR THE FIRM 

TOTCSTCJ) = TOTAL COST FOR PERIODS UP TO CJ) 

UTCSTCJ) = COST PER UNIT FOR LOT CtACERING REQUIREMENTS UP TO CJ) 
QOPT » OPTIMAL LOT SIZE FOR MINIMUM CUTCST) 

CSTMIN = MINIMUM TOTAL COST PER UNIT 

EXTIME = EXPIRY TIME IN MONTHS 

ERROR = FORECAST VARIANCE AROUND DEMAND 

ANNUAL * YEARLY DEMAND IN UNITS 

ICLASS = INDICATES THE CLASS OF THE ITEM 


DIMENSION DC4,12) f AVERGEC12),YEARAVC3),BASEC 12),RATI0C 12) 
DIMENSION FRCSTC24),T0TCSTC12),UTCSTC 12),NAMECA0) 

INPUT 0 AT A 


DATA K,N,ALPHAfRTRATEtSETUP/3f I2t0.33f 12. #40./ 
DATA T1MEAN,T2MEAN,VART1,VART2/2.#7.#1.#2./ 


DATA ISTAR/IH*/ 

KOUNT = 1 
PRINT 1000 

ICOO F0RMATC/t2Xtll5ClH-n 
PRINT 1010 

1010 FORMAT C/f2Xt*ITEM NO.*f lXt*NAME OF THE ITEM*t24X,*0PTlMAL*f 5Xf 

1 *RE0RDER*,6X,*UNIT*t4Xt*AM0UNT IN*,4X,*N0. OF*, 3X,*NEXT ORDER*,/ 

2 51X,*0RDER*,7X,*LEVEL*,8X,*PRICE*, 3X,*RS* TO BE*, 4X,*PERI00S*, 2X 

3 *T0 BE PUT*,/,51X,*QUANTITY*,17X,*IN RS.*,2X,*INVeSTED*t5X,*C3VE 
4ED*,2X,*IN MONTH*) 

PRINT 1000 

1 READ 1001,{NAMECJ),J»1,40),PRICE,EXTIME, ICLASS,NCURNT 

1001 FORMAT<40A1,F4«2,F3.0,2I3) 

IFKNAMECD.EQ.ISTARI.AND. (NAME(4D).EQ*ISTAR) ) GO TO 180 
IF ( ICLASS .EQ, 1 ) SERVCE = 0.95 

IF C IGLASS .EQ. 2 ) SERVCE = 0.99 

IF ( ICLASS .EQ. 3 ) GO TO 110 

READ 1002, T3MEAN,VART3 ‘ 

1002 F0RMAT<2F6.0) 

READ 1003, ( {D(K1, J) , J=1,N),K1=1,K) 

1003 F0RMATC12F6.0) i 

C . I 

C MONTHLY AVERAGES ARE COMPUTED IN THIS PORTION 

C ( 

00 15 J = 1,N 

SUM = 0.0 I 

DO 10 I = 1,K 

SUM = SUM + D(I,J) I 

10 CONTINUE I 

AVERGECJ) * SUM / FLOATCKl j 

15 CONTINUE 

C 

C COMPUTING THE TREND/ MONTH 

C 

DO 25 I * l,K 
SUM = 0.0 
DO 20 J = 1»N 
SUM a SUM + DCI,J) 

20 continue 

YEARAV(I) = SUM / FLOATINI 
25 CONTINUE 

TREND = (YEARAV(K)-YEARAVf 1) )/FLOAT(N*rK-ll) 

NUP * N/2 

C (N) MUST BE AN EVEN NUMBER 

C ■ : 

C DEVELOPING THE BASE SERIES AS FOLLOWS 

C ; 


DO 35 J = ItNUP 

CJ = FLOAT (J1 - 0.5 

CORFCT = CJ * TREND 

N1 = NUP - J + 1 

BASECNl) = AVERSE INI) + CORECT 

N2 = NUP + J 

BASECN2) = AVERGECN2) - CORECT 
35 CONTINUE 
C 

C COMPUTING AND SMOOTHING RATIOS AND FORECAST ERROR 

C — — — .... 

c 

RATIOCI) a 0(K»1)/BASE( 1) 

RBAR = RATIO(l) 

FRCSTf2) a BASE(2} * RBAR 

ERROR a 0.0 

DO 50 J a 2, N 

RATIO(J) a 0{K,J)/BASEC J) 

RBAR a ALPHA=«<RATI0( J) + (1.0-ALPHA) * RBAR 
NP a J + 1 

IF ( J.EQ.N ) NP a 1 
FRCST{J+1) a BASE(NP) ♦ RBAR 
ERROR a ERROR + ( C FRCST( J)-D( K, J )) /30.0)**2 
50 CONTINUE 

IF (NCURNT.EQ. 12) GO TO 200 
IF (NCURNT .EQ. 0 ) GO TO 70 
KC U ~ K + 1 
FRCST(l) a FRCSTCN+1) 

READ 1002 f (0(KCU,J),J=1, NCURNT) 

DO 60 J a 1, NCURNT 
RATIO(J) a D(KCU. J)/BASE( J) 

RBAR a ALPHA*RATI0(J) + ( 1.0-ALPHA )*RBAR 
ERROR a ERROR + ((FRCST(J) - DCKCUf J ) )/30.0) **2 
FRCST(J+l) a RBAR * BASE(J+1) 

60 CONTINUE 

ERROR a ERROR/FLO AT (N+NCURNT-l) 

C 

C DYNAMIC LOT-SIZING BY PART PERIOD ALGORITHM (MODIFIED P.P.A.) 

C 

70 CARCST a (PRICE*RTRATEI/1100.*FL0AT(N)) 

FRCST(l) a RBAR * BASEC NCURNT+1 ) 

CUANT a FRCST(I) 

TOTCST(l) a SETUP + CARCST*(0.5)*FRCST(1) 

UTCST(l) a T0TCST(1)/QUANT 
: J = 2 

75 N8 a NCURNT + J . 

NC a NB - N ■ . V:/' 

IF(NC.EQ.l) N8 a NC 

FRGST(J) a RBAR*BASE(NB) + FLOAT(J-1)*TRENO 



(-•OOOOOOOOOOOOOOOOOOO 


IIT = J-1 

TOTCSTCJ) = TOTCST(IIT) + tFLOAK J )-0.5)*FRCST{ J)*CARCST 
QUANT = QUANT + FRCSTtJI 
UTCSTCJJ = TOTCSTC4) f QUANT 
IF (UTCSTt J).GE.UTCSTtnT) ) GO TO 80 
J = J + 1 

IF CFLOATC J),GT.EXTIHE) GO TO 85 
GO TO 75 

80 QOPT = QUANT - FRCST(J) 

CSTMIN = UTCSTUIT) 

GO TO 90 
85 J = J - 1 

QOPT = QUANT 
CSTMIN = UTCSTCJ) 

90 IQOPT = IFIX( QOPT + 1. ) 

QOPT = FLOAT! IQOPT) 

QOPT = FLQATCIQOPT) 

AMOUNT = PRICE * QOPT 
ICOVER = IIT 
MONTH = NB 
GO TO 98 

THIS PORTION CALCULATES REORDER LEVEL 


TIMEAN * MEAN OF THE DISTRIBUTION OF !T1) 

T2MEAN = MEAN OF THE DISTRIBUTION OF CT2) 

T3MEAN = MEAN OF THE DISTRIBUTION OF <T3) 

VARTl = VARIANCE OF THE DISTRIBUTION OF CTl) 

VART2 = VARIANCE OF THE DISTRIBUTION OF (T2) 

VART3 = VARIANCE OF THE DISTRIBUTION OF ITS) 

TTOTAL = TOTAL LEAD TIME IN DAYS (T1+T2+T3) 

VARTOL = VARIANCE OF TOTAL LEAD TIME 

SERVCE = DESIRED SERVICE LEVEL 

DLEAD = AVERAGE DEMAND DURING THE LEAD TIME 

BUFFER = TIME IN DAYS FOR WHICH BUFFER STOCK IS TO BE KEPT 

REDDER * REORDER LEVEL 


10 IF < NCURNT .EQ. 12 ) GO TO 200 
READ 1 005 t ANNUAL, TTOTAL, BUFFER 
1005 F0RMAT<3F6«0) 

Q = SQRT((2.0*100.«ANNUAL*SETUP)/CPRICE*RTRATE)) 
IQ = IFIX(Q + 1.) 

Q = FLOAT! IQ) 

TIME = IQ+FLOATIN) 1/ANNUAL 
IF !TIME -GT. EXTIME ) GO TO 120 
QOPT = Q 

ICOVER » IFIXITIHE) 


noon no 


120 


GO TO 125 

Q = (ANNUAL^EXTIMH/FLOAKN) 

IQ = IFIX(Q > 1. Kic ; - . 

QOPT = FLOAT! IQ) " 

ICDVER = I Fixe EXTIME) 

125 NC = NCURNT + ICOVER - N 

MONTH = NCURNT IGOVER +1 
IF t NC .GE. 0 ) MONTH * NC +1 
AMOUNT a PRICE * QOPT 

REODER = IFIX!ICANNUAL*(TT0TAL+BUFFER))/365. ) + !• ) 
GO TO 130 

98 VARFOR = ERROR 

TTCTAL = TIMEAN + T2MEAN + T3MEAN 
VARTOL = VARTl + VART2 + VART3 


OMORAT = CeOPT)/! 30, ♦FLOAT! ICOVER ) ) 

DLEAO * DMDRAT ♦ TTOTAL 

VRDLED = VARFOR^!TTOTAL^^2) + V ART OL*! DMDRAT ♦♦2) + VARTOL ♦VARFOR 
REORDER LEVEL BY TCHEBYCHEFF INEQUALITY 


CONFID * SCRT! l./!4.5^!l,-SERVCE)) ) 

REDDER = IFIXCDLEAD + CONFID^SQRT! VRDLED) + 1.) 

130 CONTINUE 

PRINT 1020,K0UNT,!NAME! J), J=l,40)f QOPT, REDDER* PRICE, AMOUNTf 
1 ICCVER,MONTH 

1020 FORMAT!/,2X,I3,5X,40A1,F8.0,4X,F8.0,3X,F8.2, 2X,F8.2,7X,I2, 7X,I 2) 

PRINT 1000 
GO TO 150 

200 PRINT 2003 

2003 FORMAT!/, 20X, ♦YEAR ENDED - PLEASE UPDATE YOUR FILES^) 

GO TO 160 

180 PRINT 2005 

2005 FORMAT!//, 35X,^END OF DATA FILE^) 

GO TO 160 

150 KOUNT = KOUNT + I 

GO TO 1 

160 STOP 

END 


